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The  hypothesis  that  presence  of  the  pcriplasmic  nitrate  reductase  during  aerobic 
growth  minimizes  or  eliminates  the  diauxic  lag  was  tested.  Indeed,  little  or  no  lag  is 
observed  when  bacteria  known  to  express  the  pcriplasmic  enzyme  (Paracoccus 
/tantoirophus)  arc  subjected  to  an  anoxic  environment  after  a period  of  aerobic  growth. 
Long  lags  were  still  observed  when  aerobically  grown  cultures  of  bacteria  lacking  the 

There  was  some  indication  that  more  oxidized  substrates  led  to  longer  diauxic  lag  lengths 
for  Paracoccus  panlolrophus,  although  no  long  lags  were  observed  for  any  of  the 


The  pcriplasmic  nitrate  reductase  is  not  expressed  in  all  denitrilying  bacteria. 


the  use  of  nitrate  as  electron  acceptor.  A simple  model  has  heen  formulated  to  simulate 
diauxic  growth  of  denitrifying  bacteria  that  only  have  the  membrane-bound  enzyme.  It 
successfully  depicts  the  observed  lags  when  a pure  culture  of  Pseudomonas  denitrificuns 
switches  from  oxygen  to  nitrate  as  terminal  electron  acceptor.  The  model  includes  the 
effect  of  carbon  substrate  limitation  and  length  of  aerobic  phase;  and  docs  not  run  into 
problems  when  switching  from  anoxic  to  aerobic  conditions,  unlike  earlier  models  of 
diauxic  growth.  In  order  to  lit  growth  data  of  bacteria  that  also  have  the  pcriplosntic 
nitrate  reductase,  a preliminary  model  was  devised  that  includes  the  contribution  of  both 


CHAPTER  1 
INTRODUCTION 

Removal  of  nitrogen  compounds  from  wastewater  is  very  important  to  satisfy 
increasingly  stringent  effluent  quality  requirements.  There  arc  several  reasons  why  this  is 
a necessary  process.  Presence  of  nitrogen  contributes  to  eutrophication  of  bodies  of 
water,  exerts  an  oxygen  demand  that  can  cause  harmful  effects  on  aquatic  life  and  the 
environment,  and  is  toxic  in  some  forms  (Ramalho,  1983).  Biological  processes  are  the 
preferred  means  of  removing  nitrogen  from  wastewater  (Grady  et  al.,  1999). 

Biological  nitrogen  removal  in  wastewater  treatment  plants  is  achieved  through 
bacterially  mediated  oxidation  of  ammonia  to  nitrate;  and  reduction  of  nitrate  to  nitrite, 
nitrous  oxide,  nitric  oxide,  or  dinitrogen.  Accumulation  of  any  of  the  intermediates  (i.c., 
nitrite,  nitric  oxide,  or  nitrous  oxide)  is  undesirable.  Biological  oxidation  of  ammonia  to 
nitrate  is  termed  nitrification  and  is  carried  out  by  aerobic,  autotrophic  bacteria  (e.g„ 
Nitrosomonas  and  Nitrohacler).  Biological  reduction  of  nitrate  to  dinilrogen  is  known  as 
complete  denitrification.  Denitrification  is  carried  out  by  hctcrotrophic,  facultative 
anaerobic  bacteria  (e.g..  Pseudomonas  denunficans.  Paracoccus pamotrophus. 
Pseudomonas  stutzeri ).  Although  denitrification  under  aerobic  conditions  hus  been 

nitrate  present)  conditions  (Ticdjc,  1988).  Therefore,  to  achieve  nitrogen  removal  in  a 
single-sludge  system,  a mixed  population  of  nitrifying  and  denitrifying  bacteria  must  be 
exposed  to  alternating  aerobic  and  anoxic  conditions  (Ramalho.  1983).  In  order  for  this 
process  to  be  efficient,  the  bacteria  must  be  able  to  adapt  quickly  to  these  changes. 


Several  investigators  (Gouw  ct  al.,  200 1 ; Liu  et  al.,  1 998a.  1 998b;  Lisbon  ct  al., 
2002)  have  shown  that  when  bacteria  change  from  dissolved  oxygen  to  nitrate  as  electron 
acceptors,  they  can  experience  a period  of  little  or  no  growth.  Such  a period,  when  it 
occurs  between  two  periods  of  exponential  growth,  was  termed  diauxic  lag  by  Monod 
(1942).  The  same  terminology  is  used  in  this  study.  Lag  periods,  since  they  reduce  the 
overall  rate  of  nitrogen  removal,  are  undesirable  in  nitrogen-removing  wastewater 
treatment  processes.  It  is  believed  that  during  this  lag  bacteria  will  produce  the 
membrane-bound  nitrate  reductase,  an  enzyme  that  is  inhibited  by  oxygen  and  is 

Some  researchers  have  reported  simultaneous  utilization  of  nitrate  and  oxygen 
(Robertson  and  Kuencn,  1 984)  by  Thiosphaera  paniotropha  (now  known  as  Paracoccus 
pantolrophus).  These  bacteria  arc  known  to  have  a pcriplasmic  nitrate  reductase  (Nap) 
that  is  not  inhibited  by  oxygen.  It  has  been  suggested  that  a possible  role  for  this  enzyme 
is  to  aid  in  the  transition  from  environments  where  oxygen  is  present  to  environments 
where  it  is  absent  (Moreno-Vivita  et  al.,  1 999).  If  this  is  the  case,  there  is  a possibility 

been  suggested  that  Nap  is  used  for  redox  balancing,  and  that  the  activity  ofNap 
increases  with  the  extent  of  reduction  of  the  carbon  substrate  used  for  growth 

containing  bacteria  will  decrease  as  the  reduction  state  of  the  carbon  substrate  increases. 


CHAPTER  2 

REVIEW  OF  LITERATURE 

Nitrogen  Removal 

Nitrogen  removal  can  be  achieved  by  several  physico-chemical  as  well  as 
biological  processes.  Reactions  generating  nitrogen-containing  gases  catalyzed  by 
chemical  agents  are  known  as  chcmodcnitrificalion.  An  example  of  these  is  the  acid- 
catalyzed  destruction  of  nitrite,  in  which  the  predominant  product  is  usually  NO.  which  is 
undesirable  (Ticdje,  1 988).  The  gas.  NO.  is  a very  minor  product  of  biological 
denitrification,  which  is  one  of  the  reasons  why  biological  denitrification  is  a very 
attractive  method  for  nitrogen  removal.  Furthermore,  physico-chemical  technologies  are 
expensive  and  difficult  to  maintain  (Grady  et  al.,  1999). 

Biological  nitrogen  removal  from  municipal  wastewater  requires  both  nitrification 
and  denitrification.  In  nitrification,  ammonia  is  oxidized  to  nitrite  (NOj")  in  the  presence 
of  oxygen  by  specific  bacteria  ( Nitrosomonas ) and.  in  turn,  nitrite  is  oxidized  to  nitrate 
(NOf)  in  the  presence  of  oxygen  by  other  specific  bacteria  ( Nitrohacter ).  In 
denitrification,  nitrate  is  converted  to  nitrogen  gases  (NO.  NjO.  N;)  in  the  absence  of 
molecular  oxygen  (Ramalho,  1983). 

Many  of  the  bacteria  that  can  denitrify  are  facultative  anaerobes,  such  as 
Pseudomonas  denilriflcans.  These  bacteria  can  grow  under  both  anoxic  and  aerobic 
conditions  (Komaros  and  Lyberatos.  1998).  However,  oxygen  respiration  is  more 
energetically  efficient  than  denitrification  (respiration  of  nitrogen  oxides),  so  facultative 


bacteria  usually  utilize  nitrogen  oxides  only  in  the  absence  of  oxygen  (Oh  and 
Silverstein,  1999b). 

A number  of  approaches  were  taken  when  developing  the  basis  for  biological 
nitrogen  removal.  One  involved  the  use  of  a series  of  separate  suspended  growth  systems 
to  accomplish  removal  of  organic  matter  and  nitrogen  in  a step-wise  fashion.  Organic 
matter  was  removed  in  the  first  step,  nitrification  occurred  in  the  second,  and 
dennnficauon  occurred  in  the  third  (with  addition  of  an  external  carbon  source).  These 
systems  involved  high  construction  and  operating  costs.  Another  approach,  involving  the 
incorporation  of  aerobic  zones  for  nitrification  and  anoxic  zones  for  denitrification  in  a 
single  system,  with  carbon  oxidation  occurring  in  both  zones,  was  developed.  It  was 
referred  to  as  single-sludge  nitrogen  removal  (Grady  el  al„  1 999).  In  these  systems,  the 
bacterial  population  must  be  alternately  exposed  to  aerobic  and  anoxic  environments. 
This  forces  the  biomass  to  switch  between  electron  acceptors,  using  oxygen  during  the 
aerobic  phase  and  nitrate  in  the  absence  of  oxygen. 

The  required  sizes  of  facilities  for  biological  denitrification  depend  on  the  kinetics 
of  denitrifying  bacteria.  Several  factors  may  affect  denitrification  kinetics.  For  example, 
the  quantities  of  nitrogen  oxide  reductases  (enzymes  required  for  utilization  of  nitrogen 
oxides)  are  affected  by  induction  (if  present)  by  NO,'.  NO2',  NO.  and  N2O,  or  repression 
by  O2.  O2  might  decrease  enzyme  activity  due  to  competition  for  electrons  or  by  enzyme 
inactivation.  Limited  electron  donor  (shortage  of  available  carbon)  might  also  reduce 
enzyme  activity;  as  might  other  factors,  such  as,  pH,  temperature,  and/or  toxic  chemicals 
(Tiedje,  1988). 


Regulation  of  Enzymes  Involved  in  Denitrification 


Reduction  of  nitrate  in  anaerobic  environments  occurs  mostly  by  two 


dissimilatory  processes:  respiratory  denitrification  and  dissimilatoiy  nitrate  reduction  to 
ammonium.  Although  assimilatory  nitrate  reduction  could  also  occur,  presence  of  high 
concentrations  of  ammonium  and  organic  nitrogen,  which  commonly  exist  in  anaerobic 
environments,  repress  this  process.  Dissimilatoiy  processes  differ  from  assimilatory 


processes  by  the  fact  that  the  reduced  nitrogen  is  not  used  by  the  cell.  Nitrogen  oxides  in 
dissimilatory  processes  serve  as  electron  acceptors  for  the  cell's  metabolism  (Ticdje. 
1988). 


Oxygen  regulates  denitrification  by  repressing  enzyme  synthesis  or  by  inhibiting 
enzyme  activity  (Tiedje.  1 988).  It  is  generally  accepted  that  the  respiratory  nitrate 
reductase  and  other  respiratory  nitrogen  oxide  reductases  are  repressed  or  inhibited  by  Oj 
to  different  extents  (Tiedje.  1988;  Payne.  1981).  The  later  reductases  (for  NO  and  NjO) 
in  the  denitrification  sequence  are  somewhat  more  Or  seasitive  than  the  earlier  reductases 
(Knowles.  1982).  Bonin  et  al.  (1989)  studied  the  effects  of  oxygen  on  each  step  of 
denitrification  by  Pseudomonas  nautica,  and  observed  that  every  enzyme  associated  with 
denitnficauon  is  affected  differently  widt  respect  to  oxygen  concentration.  Nitrate 
reductase  turned  out  to  be  less  sensitive  to  oxygen  than  the  other  nitrous  oxide  reductases 
were.  At  oxygen  concentrations  greater  than  4.05  mg/L,  nitrate  reductase  was  completely 
inactivated,  while  nitrite  and  nitrous  oxide  reductase  activities  were  blocked  at  2.15  mg/L 
and  0.25  mg/L  respectively.  McKcnncy  cl  al.  ( 1 994)  studied  the  effects  of  oxygen  on  the 
kinetics  of  denitrification  by  Pseudomonas fluorescent  Their  experiments  showed  dial 
oxygen  inhibited  all  steps  in  the  denitrification  process.  The  inhibition  was  rapid  afier 
imroducing  oxygen  and  was  reversible  after  restoring  anoxic  conditions.  The  transient 


effects  suggested  that  the  reduction  of  NO;'.  NO  and  N2O  arc  each  inhibited  to  increasing 
extents  by  oxygen.  However,  completely  anaerobic  conditions  were  not  necessary  for 
denitrification  to  occur. 

The  reduction  of  nitrate  to  nitrite  is  catalyzed  by  three  different  bacterial  enzymes: 
the  cytoplasmic  assimilatory  (Nas),  the  membrane-bound  respiratory  (Nar),  and  the 
pcriplasmic  dissimilatory  (Nop)  nitrate  reductases. 

Regulation  of  Nas  by  ammonium  takes  place  in  two  ways:  repression  by  the 
general  nitrogen  regulatory  system  (Ntr).  which  regulates  synthesis  of  most  enzymes 
required  for  use  of  alternate  nitrogen  sources;  and  inhibition  of  nitrate  transport,  avoiding 
induction  of  nitrate  rcduclase  t Moreno- Viviun  ct  al..  1999). 

A membrane-bound  nitrate  rcduclase  (Nar)  has  been  identified  in  several  nitrate- 
utilizing  bacteria.  The  activity  of  Nar  in  bacteria  is  inhibited  by  oxygen.  Observations 
with  Paracoccus  denitrificans  indicated  that  this  is  not  a direct  inhibition.  Instead,  it  is 
exerted  via  the  redox  state  of  a cellular  component  controlling  the  movement  of  nitrate 
across  the  cytoplasmic  membrane  to  the  active  site  of  the  reductase  ( Alcfoundcr  and 
Ferguson.  1980).  Nar  is  predominantly  expressed  during  growth  under  anaerobic 
conditions  with  nitrate  present  and  is  only  operative  in  the  absence  of  oxygen  (Bell  et  al.. 
1990).  However,  Krul  and  Vecningen  (1977)  (bund  that  in  Alcaligenes  strain  15,  Nor 
only  partially  lost  activity  in  the  presence  of  oxygen.  Thus,  it  is  not  conclusive  that  in  all 
microorganisms  Nar  is  inoperative  under  aerobic  conditions.  Azide,  cyanide,  and 
thiocyanate  have  been  found  to  inhibit  this  enzyme  (Hochstcin.  and  Tomlinson.  1988). 

The  second  dissimilatory  nitrate  reductase  (Nap),  which  is  located  in  the 
periplasm,  has  been  found  in  aerobically  denitrifying  bacteria.  Nap  is  predominantly 


absence  of  oxygen.  Nap  is  not  affected  by  presence  of  ammonium.  Palureau  et  al. 
(1994).  in  Ihoir  studies  of  denitrification  by  Comamonas  sp.  under  aerobic  conditions, 

of  previously  synthesized  Nap. 

Several  authors  have  suggested  that  Nap  is  a dissimilatory  enzyme  used  to  shunt 
excess  electrons  to  nitrate  for  redox  balancing  (Richardson  and  Ferguson,  1 992;  Scars  ct 
al.,  2000).  Redox  balancing  can  be  necessary  for  optimal  bacterial  growth  under  some 

substrates  in  aerobic  hetcrotrophs  (Moreno-Vivian  et  al.,  1999).  The  level  of  Nap  activity 
during  aerobic  growth  has  been  shown  to  increase  with  the  extent  of  reduction  of  the 
carbon  substrate  in  cultures  of  Thiosphaera  pantotropha.  A possible  reason  for  this  is 

regulation  leading  to  the  activity  change  in  response  to  redox  state  is  still  not  clear. 

of  two  promoters  upstream  of  the  nap  gene  cluster.  One  promoter  is  mote  active  during 
aerobic  growth  on  reduced  carbon  substrates  and  the  other  during  aerobic  growth  on 
oxidized  carbon  substrates.  The  activity  of  both  promoters  is  repressed  under  anaerobic 
conditions,  when  nar  genes  are  expressed. 

Because  Nap  is  located  in  the  periplasm,  nitrate  transport  limitation  by  oxygen  (as 
in  the  case  of  the  Nar  system)  is  not  expected.  Aerobic  denitrification  is  then 


ol.  ( 1 999)  proposed  that  Nap  is  required  to  catalyze  the  find  step  of  denitrification  in 
Pseudomonas  sp.  Strain  G-179.  They  saw  no  growth  of  a Nap-deficient  mutant  under 
denitrifying  (anoxic)  conditions,  it  is  not  clear  under  what  conditions  cultures  were 
grown  prior  to  the  anoxic  conditions.  Perhaps  there  was  simply  a lag  and  they  did  not 
wait  long  enough  for  the  cultures  to  build  up  Nar  and  begin  anoxic  growth. 

It  has  also  been  suggested  that  Nap  may  be  utilized  in  the  adaptation  to 
metabolism  in  the  absence  of  oxygen  after  transition  from  aerobic  conditions  (Moreno- 
Vivian  el  al.,  1 999).  This  latter  suggestion  is  especially  pertinent  to  wastewater  treatment 
plants  that  alternate  the  environment  of  biological  populations  between  aerobic  and 

Wamcckc-Ebcrz  and  Friedrich  (1993)  studied  the  activities  of  three  nitrate 
reductases  (Nas,  Nar.  and  Nap)  in  Alcaligenes  eutrophus  strain  H 16.  In  agreement  with 
other  studies,  they  found  that  oxygen  shuts  down  synthesis  of  Nar  and  that  ammonium 
shuts  down  synthesis  of  Nas.  Interestingly,  the  activity  of  Nap  in  this  organism  was  only 
observed  during  stationary  phase;  whereas  Richardson  and  Ferguson  (1992)  found  that,  in 
Thiosphaem panlolropha,  the  pcriplasmic  dissimilatory  nitrate  reductase  was  expressed 
during  exponential  growth. 

The  dissimilatory  nitrate  reductase  of  most  microorganisms  is  induced  by  nitrate 
under  anaerobic  conditions.  However,  sometimes  nitrate  docs  not  need  to  be  present  to 
induce  this  enzyme  (Schulpand  Stouthamcr,  1970;  Krul  and  Vceningcn,  1977). 

It  is  commonly  believed  that  oxygen  represses  the  synthesis  of  the  membrane- 
bound  dissimilatory  nitrate  reductase  (c.g.,  Wamccke-Eberz  and  Friedrich,  1 993; 
Moreno-Vivian  et  al.,  1999).  Krul  and  Vecningcn  (1977)  studied  the  synthesis  of  nitrate 


rcducuisc  under  aerobic  conditions  by  various  denitrifying  bacteria  isolated  from 
activated  sludge  and  drinking  water,  including  Alcaligenes  strains  1 5 and  N4.  Although 
chlorate  can  be  reduced  by  the  dissimilatory  nitrate  reductase  (Pichinoty.  1 965).  it  was 
shown  to  repress  the  synthesis  of  the  dissimilatory  nitrate  reductase  in  strain  1 5.  The 
nitrate  reductase  activity  measured  with  cells  prccultured  in  the  presence  of  high  amounts 
of  chlorate  is  comparable  to  the  activity  measured  in  cells  grown  at  high  dissolved 
oxygen  concentrations  without  chlorate.  This  can  be  explained  by  the  hypothesis  of 
Wimpenny  and  Cole  (1967)  that  redox  potential  is  the  regulating  factor  for  synthesis  of 
dissimilatory  nitrate  reductase.  Thus,  the  redox  potential-increasing  effect  of  chlorate 
might  be  responsible  for  the  repression.  However,  synthesis  of  dissimilatory  nitrate 
reductase  was  not  seriously  repressed  by  chlorate  in  strain  N4;  suggesting  either  a 
different  regulation  pattern,  or  the  possibility  that  the  dissimilatory  nitrate  reductase  in 
N4  is  constitutive.  It  seems  that,  at  least  for  some  bacteria,  it  is  redox  potential,  as 
opposed  to  oxygen  concentration  per  sc,  that  represses  biosynthesis  of  membrane-bound 
dissimilatory  nitrate  reductase. 


Effect  of  Oxygen  on  Nitrogen  Removal  through  Denitrification 

As  the  previous  discussion  on  enzyme  regulation  makes  clear,  oxygen  serves  as 
an  inhibitor  of  activity  or  as  a repressor  of  synthesis  of  denitrification  enzymes.  The 
existence  of  a threshold  oxygen  concentration  above  which  denitrification  cannot  occur 
has  been  suggested  because  of  the  observed  effect  of  oxygen  on  nitrogen  removal 
through  denitrification.  This  implies  a critical  concentration  instead  of  a gradual 
inhibitory  response.  Different  values  have  been  reported  and  show  that  active 
denitrifying  microorganisms  in  marine,  freshwater,  and  soil  environments  have  very  low 


thresholds  for  oxygen.  However,  wastewater  environments  that  provide  an  excess  of 


electron  donor  and  a fluctuating  status  of  nitrate  and  oxygen  might  select  for  organisms 
that  could  make  simultaneous  use  of  both  nitrate  and  oxygen  as  electron  acceptors 
(Ticdje,  1988).  This  phenomenon  has  been  termed  aerobic  denitrification. 

Aerobic  denitrification  has  been  observed  by  several  authors.  Robertson  and 
Kucnen  found  that  Thiosphaera  panlolropha  (LMD  82.5)  is  capable  of  simultaneously 
utilizing  nitrate  and  oxygen  as  terminal  electron  acceptors  in  respiration.  This  behavior 
has  been  found  in  cultures  maintained  at  dissolved  oxygen  concentrations  as  high  as  90% 
of  air  saturation  (Robertson  and  Kuencn,  1 984).  When  using  both  terminal  electron 
acceptors  simultaneously,  batch  cultures  of  Thiosphaera  panlolropha  grew  more  rapidly 
than  when  supplied  with  cither  electron  acceptor  alone.  In  two-member  chemostat 
cultures  containing  Thiosphaera  panlolropha  and  Paracoccus  Jenilrijicans  (LMD  22.21), 
a strictly  anaerobic  denilrifier.  the  first  was  relatively  unaffected  by  acrobic/anaerobic 
switches,  as  long  as  nitrate  was  present.  However,  the  latter  washed  out  when  suddenly 
switched  from  aerobic  to  anaerobic  conditions,  possibly  because  it  was  unable  to 
synthesize  its  denitrification  enzymes  quickly  enough  (Robertson  und  Kucnen.  1990). 

Strain  SGLY2  of  Comamonas  sp.  also  is  capable  of  denitrifying  under  aerobic 
conditions  (Patureau  et  al.,  1994).  According  to  their  study,  oxygen  was  used 
simultaneously  with  nitrate  for  respiration,  giving  a growth  rate  ten  times  greater  than 
that  obtained  when  only  nitrate  was  used  as  an  electron  acceptor.  Nitrogen  gas 
production  was  observed  during  aerobic  nitrate  reduction. 

Denitrification  at  positive  bulk  dissolved  oxygen  concentrations  has  been  widely 
observed  in  activated  sludge  (e.g.,  Krul.  1976;  O'Neill  and  Horan.  1995;  Munch  et  al.. 
1996;  Goronszy,  1997).  The  typical  observations,  as  reiterated  by  Oh  and  Silverstcin 


(1999a),  are  that  low  dissolved  oxygen  concentrations  decrease  the  specific 
denitrification  rale  and  this  effect  is  amplified  with  increasing  dissolved  oxygen 


concentrations.  Furthermore,  high  dissolved  oxygen  concentrations  produce  almost 


complete  inhibition  of  denitrification.  At  dissolved  oxygen  concentrations  as  high  as  2 
mg/L.  they  found  that  denitrification  continued,  although  at  a greatly  reduced  rate. 

Several  possible  explanations  exist  for  this  behavior.  One  is  the  presence  of 
aerobically  denitrifying  bacteria,  such  as  Thiosphaera panlolrophu  and  Comamonas  sp. 
strain  SGLY2,  as  described  above.  It  could  also  be  that,  even  if  oxygen  inhibits  the 
denitrifying  activity  of  bacteria,  anoxic  regions  can  exist  in  an  activated  sludge  floe  due 
to  limitation  of  oxygen  diffusion  within  the  Hoc.  These  micro-niches  allow  nitrification 
to  occur  in  the  outer  region  in  contact  with  bulk  water  DO  (dissolved  oxygen),  and 
denitrification  in  an  inner  anoxic  zone. 

Krul  ( 1 976)  did  a study  to  compare  the  oxygen  uptake  and  dissimilatory  nitrate 
reduction  by  cells  in  suspension  and  cells  in  floes.  A number  of  bacteria  were  isolated 


from  aerobically  grown  and  denitrifying  activated  sludge:  and  were  tested  for  their 
capacity  to  form  floes,  sensitivity  of  the  floes  to  dispersion,  and  ability  to  denitrify.  Floes 
of  anaerobically  grown  cells  of  Alcaligenes  strain  1 5 showed  greatly  increased  nitrate 
reduction  below  1 ,5  mg/L  oxygen  as  opposed  to  cells  in  suspension,  which  did  not  show 
increased  reduction  until  the  concentration  of  oxygen  was  close  to  zero.  This  difference 
in  behavior  become  more  evident  when  floes  and  suspended  cells  had  been  pre-aeraled  in 
a nitrogen-free  medium  for  24  hours  before  being  tested. 

Thus,  the  anaerobic,  denitrifying  inner  part  of  floes  may  be  responsible  for  the 
loss  of  nitrogen  in  activated  sludge  plants  under  aerobic  conditions.  This  explanation  has 
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oxygen  concentration  increased.  The  fully  aerobic  sample  did  not  release  nitrous  oxide. 
When  switching  from  aerobic  conditions  back  to  anaerobic  conditions,  the  activity  of 
nitrite  reductase  was  recovered.  An  experiment  was  performed  with  two  cycles  of 
anaerobic/aerobic  switches.  Although  inactivation  of  nitrite  reductase  occurred  in  both 
cycles,  the  production  of  nitrous  oxide  was  greater  in  the  second  cycle.  This  suggests 
that  inactivation  by  oxygen  of  nitrite  reductase  synthesized  during  anaerobic  prcculture  is 
reversible.  Activity  of  the  enzyme  is  lost  under  aerobic  conditions,  and  is  regained  upon 
switching  to  anaerobic  conditions.  To  prove  that  the  nitrite  reductase  from  the  initial 
anaerobic  phase  was  responsible  for  the  nitrile  reducing  activity  in  the  second  anaerobic 
phase,  chloramphenicol  was  added  just  before  shifting  from  aerobic  to  anaerobic 
conditions.  This  would  Irak  all  protein  synthesis  in  the  ensuing  anaerobic  period.  The 
evolution  of  nitrous  oxide  in  the  second  anaerobic  phase  was  not  affected.  The 
persistence  of  nitrite  reductase  was  examined  for  various  durations  of  aerobic  periods. 
The  remaining  activity  (measured  under  anaerobic  conditions)  was  practically  the  same 
after  up  to  30  hours  of  exposure  to  aerobic  conditions  (Kawakami  et  al.,  1985). 

Otte  ct  al.  (1996)  studied  denitriliftcation  (nitrite  reduction)  by  AlcaliRenes 
faecalh  in  continuous  reactors  operating  at  different  oxygen  concentrations  or  switching 
between  aerobic  and  anaerobic  conditions.  Nitrous  oxide  (N2O)  was  produced  under  all 
conditions,  from  fully  anaerobic  to  fully  aerobic.  Nitrite  reductase  activities  under  fully 
anaerobic  conditions  were  much  higher  than  under  low  dissolved  oxygen  conditions, 
indicating  higher  nitrite  consumption  under  fully  anaerobic  conditions.  Yet.  the  highest 
relative  N.O  emissions  were  observed  under  low  dissolved  oxygen  concentrations. 

During  the  transition  front  anaerobic  to  aerobic  conditions.  Nd)  immediately  started  to 


after  a high  (1 30  mg/L  as  N)  i 


Baumann  el  al.  (1997)  studied  the  effect  of  repeated  alternating  acrobic/anoxic 


growth  conditions  on  the  development  of  stable  denitrification  activity  (i.e.,  nitrogen  gas 
was  the  only  product  released)  with  Paracoccus  denilriflcans  and  with  activated  sludge  in 
continuous  reactors.  They  found  that  the  length  of  the  anaerobic  period  influenced 
biosynthesis  of  the  denitrification  enzymes.  Long  anaerobic  periods  led  to  stable 
denitrification  activity  producing  only  dinitrogen.  Short  anaerobic  periods  led  to  nitrite 
accumulation  in  the  medium.  It  was  concluded  that  the  different  denitrification  enzymes 
in  these  cultures  arc  not  induced  simultaneously,  but  sequentially.  They  found  similar 
results  when  studying  the  effects  on  activated  sludge.  Additionally,  they  found  that  once 
the  denitrification  enzymes  are  synthesized,  they  are  not  degraded  during  aerobic  growth 
conditions.  Their  catalytic  activity  is  inhibited. 

The  effect  of  alternating  aerobic  and  anoxic  conditions  on  denitrification  by 
activated  sludge  was  also  investigated  by  Oh  and  Silverstein  (1999b).  After  several  hours 
of  anoxic  denitrification,  mixed  liquor  was  aerated  for  variable  amounts  of  time  and  then 
returned  to  anoxic  conditions.  The  specific  denitrification  rate  in  the  second  anoxic  phase 
was  decreased  in  every  cose.  In  order  to  separate  the  effects  of  dissolved  oxygen  alone 
from  the  effects  of  aerobic  growth  on  denitrification,  these  investigators  performed 
experiments  where  substrate  was  not  added  during  the  aerobic  phases.  Aeration  without 
carbon  for  long  periods  of  time  had  no  effect  on  the  anoxic  specific  denitrification  rate. 

As  mentioned  previously,  oxygen  has  two  possible  modes  of  action  in  regulating 
denitrification:  repression  of  synthesis  of  nitrogen  oxide  reductases  and  inhibition  of 
reductase  activity.  Significant  lags  when  switching  from  aeration  to  anoxic 
denitrification  are  likely  to  be  associated  with  repression,  whereas  inhibition  of  enzyme 


hum,  synthesis  of  fl-gala 


longer  in  (.'fleet,  allowing  synthesis  of  p-galactosidasc.  A Per  a while  (when  sufficient  |i- 
galactosidase  has  been  synthesized)  growth  on  lactose  begins.  The  time  required  for 
accumulation  of  P-galactosidasc  corresponds  to  the  diauxic  lag. 

Enzymes  subject  to  catabolitc  repression  require  binding  of  the  catabolitc 
activator  protein  (CAP)  to  DNA,  before  RNA  polymerase  can  bind  and  transcription  can 
take  place.  The  CAP  protein  binds  to  DNA  only  alicr  binding  cyclic  adenosine 
monophosphate  (cyclic  AMP).  Glucose  inhibits  the  synthesis  of  cyclic  AMP  and 

binding  of  RNA  polymerase  to  the  promoter. 

mixtures  of  substrates  resulting  in  a diauxic  growth  pattern  is  most  commonly  observed. 
However,  the  simultaneous  use  of  different  carbon  sources  lias  also  been  found.  “The 
preferential  growth  pattern  (sequential  utilization)  occurs  because  synthesis  of  the 
enzymes  that  catalyze  the  transport/catabolism  of  the  less  preferred  substrate  is 
completely  suppressed  in  the  presence  of  the  preferred  substrate.  On  the  other  hand, 
simultaneous  utilization  of  the  substrates  occurs  when  suppression  of  enzyme  synthesis  is 
incomplete.  Thus,  the  preferential  growth  pattern  occurs  because  the  enzyme  catalyzing 
the  uptake  of  the  less  preferred  substrate  becomes  extinct;  the  simultaneous  growth 
pattern  occurs  because  the  enzymes  for  both  substrates  coexist."  (Narang,  1 998).  An 
analogy  can  be  drawn  between  the  sequential  versus  simultaneous  utilization  of  electron 
donots  (carbon  substrates)  and  the  sequential  vetsus  simultaneous  utilization  of  electron 


acceptors  (e.g.,  oxygen,  nitrate,  nitrite).  (Reminder:  Aerobic  denitrification  is  the 
simultaneous  utilization  of  oxygen  and  nitrate.) 

Biphasic  growth  due  to  a change  in  electron  acceptors  (nitrate  to  nitrite)  was  first 
studied  by  Kodama  ct  al.  (1969).  They  observed  that  growth  of  Pseudomonas  stutseri  in 
the  presence  of  20  mM  NOj'  (280  mg/L  os  N)  or  higher  (with  oxygen  absent)  results  in 
two  exponential  phases  separated  by  an  intermediate  phase  in  which  the  growth  is  very 
low.  In  the  first  exponential  phase,  nitrate  is  reduced  to  nitrite;  and  in  the  second,  nitrile 
is  reduced  to  nitrogen  gas.  They  explain  that  'The  mechanism  behind  this  behavior  is  one 
in  which  nitrate  inhibits  nitrite  reduction  in  two  ways:  by  repressing  formation  of  the 
nitrite-reducing  system  and  by  competing  with  nitrite  as  an  electron  acceptor  in  the 
oxidative  function  of  cells"  (Kodama  et  al.,  1969,  p.  864). 

Another  observation  of  diauxic  lag  resulting  from  a transition  from  nitrate  to 
nitrite  as  electron  acceptor  was  reported  by  Komaros  ct  al.  (1996).  They  conducted  batch 
experiments  in  which  Pseudomonas  denilrificans  were  first  grown  in  the  presence  of 
nitrite,  with  oxygen  and  nitrate  absent.  Nitrate  was  then  injected,  to  give  a medium 
concentration  of  33  mg/L  as  N.  Immediately  after  the  injection,  nitrite  started 
accumulating;  and  after  nitrate  was  exhausted,  a short  lag  phase  was  observed,  after 
which  growth  resumed  on  nitrite. 

The  first  observation  of  a diauxic  lag  experienced  by  bacteria  switching  from 
oxygen  to  nitrate  as  a terminal  electron  acceptor,  was  reported  by  Waki  et  al.  (1980). 
Although  they  did  not  use  the  terms  biphasic  growth  or  diauxic  lag  when  discussing  their 
results,  the  bacterial  cultures  studied  clearly  exhibited  this  phenomenon. 


r (Liu  et  al. 


denitrificans  (ATCC  13867),  denitrifying  bacteria  that  were  originally  isolated 

1998b).  The  length  of  the  aerobic  growth  phase  was  found  to  significantly  affect  the 
length  of  the  diauxic  lag.  Longer  lags  resulted  from  longer  aerobic  growth  phases. 
Contrary  to  the  rapid  reduction  of  nitrate  observed  by  Waki  et  al.  (1980)  immediately 
after  Paracoccus  denitrificans  switched  ftom  aerobic  to  anoxic  conditions,  nitrate 

increased  rapidly  as  exponential  growth  began  to  take  place.  Presence  or  absence  of 
nitrate  during  the  aerobic  phase  had  a marginal  effect  on  lag  length  (Liu  et  al.  1998b). 
The  influence  of  preculture  conditions  on  the  diauxic  lag  of  Pseudomonas 

studied  by  Gouw  et  al.  (2001).  (Here,  preculture  refers  to  a phase  of  growth  before  the 
start  of  sequential  aerobic  and  anoxic  growth  phases).  When  bacteria  were  precullured 

were  precultured  in  the  presence  of  both  oxygen  and  nitrate.  Presence  of  nitrate  during 
long  aerobic  phases  somewhat  compensated  for  the  absence  of  nitrate  during  precultute 
phase.  In  experiments  with  anoxic  ptecullure,  short  logs  were  usually  observed  when 
nitrate  was  absent  from  the  aerobic  phase.  Even  shorter  lags,  or  no  lags  at  all,  were 
observed  when  nitrate  was  present  during  the  aerobic  phase. 


i and  Lybcratos  (1998)  found  i 


authors,  it  became  known  that  the  conditions  in  their  experiments  were  similar  to  those 
under  which  Gouw  et  ul.  (2001 ) found  no  lags  upon  the  switch. 

Lisbon  et  al.  (2002)  studied  the  effect  of  dissolved  oxygen  concentration  in  the 
aerobic  phase  on  the  ensuing  lag  of  Pseudomonas  denitrificans  when  switching  to  anoxic 
conditions.  Exposure  of  the  bacteria  to  low  dissolved  oxygen  concentrations  during  the 
aerobic  phase  resulted  in  significantly  shorter  lags  than  did  exposure  to  higher 
concentrations  approaching  air  saturation.  The  cultures  grown  under  low  dissolved 
oxygen  concentrations  generally  grew  faster  during  the  subsequent  anoxic  phase. 

Casasus-Zambrana  (2002)  studied  the  effects  of  dissolved  oxygen  and  aerobic  growth 
on  the  ensuing  diauxic  lag  and  anoxic  specific  growth  rate  of  Psemlomonas  denilrificans. 
A lag  was  more  often  observed  when  the  bacteria  are  subjected  to  an  anoxic  environment 
after  a period  of  aerobic  growth  than  after  a period  of  exposure  to  air  without  growth. 
Presence  or  absence  of  substrate  during  the  aerobic  phase  had  no  effect  on  the  specific 
growth  rate  in  the  ensuing  anoxic  phase. 

In  order  to  optimize  the  design  of  wastewater  treatment  facilities,  reliable  models 
must  be  formulated.  Successful  models  using  a cybernetic  approach  have  been 
developed  to  simulate  the  diauxic  lag  experienced  by  bacteria  switching  between  electron 
donors  (e.g,,  Ramkrishna,  1982;  Ramkrishna  etal.,  1984;  Kompalaetal.,  1986).  These 
models  propose  that  bacteria  regulate  enzymes  in  a way  that  will  maximize  their 
instantaneous  specific  growth  rate.  Narang  et  al.  ( 1 997a)  studied  whether  or  not  these 
models  could  predict  behaviors  in  systems  with  simultaneous  utilization  of  substrates,  in 


tic  lag  lakes  place.  They  found  lhat  Ihe  model  predicts  pn 


utilization  patterns  for  all  initial  conditions. 

Traditional  models  of  single-sludge  processes,  such  as  the  International  Water 
Association  (I  WA)  Activated  Sludge  Models  (ASM)  1 . 2, 2d.  and  3 (Henze  ct  al.,  2000) 

used  the  cybernetic  approach  ofKompalactal.  (1986)  to  modify  ASM  No.  1 so  that  it 

conditions.  The  model  they  obtained  was  later  modilied  to  show  the  considerably  longer 
lags  observed  when  Pseudomonas  denitrificans  switched  between  oxygen  and  nitrate  us 
electron  acceptors,  as  well  as  the  dependence  of  this  lag  on  the  length  of  the  aerobic 
growth  phase  (Liu  et  al..  1998b).  However,  it  runs  into  trouble  when  growth  switches 
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This  variable  is  eliminated,  since  the  presence  or  absence  of  nitrate  docs  not  affect  the 
synthesis  of  oxygenase. 

The  biomass  decay  rate  is  given  as  in  the  IWA  Activated  Sludge  Models  by 
Equation  4-5. 


Here,  b is  the  biomass  specific  decay  rate.  The  rale  of  decay  of  reductase  is  given  by 
Equation  4-6. 


Here  pm  is  the  reductase  specific  decay  rate.  Analogously,  the  rate  of  decay  of 
oxygenase  is  given  by 


(4-7) 


where  Po  is  the  oxygenase  specific  decay  rale. 

For  non-limiting  nitrate  and  substrate  concentration  and  zero  dissolved  oxygen 
concentration,  the  muximum  value  of  em  results  in 


which  yields  Equation  4-9. 


(4-9) 


Analogously,  the  maximum  value  of  e„  (eojma)  under  optimal  conditions  (non-limiting 
dissolved  oxygen  and  substrate  concentration)  is  Equation  4-10. 


(4-10) 
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Experimental  Methods 

All  experimental  methods  were  performed  under  aseptic  conditions.  Freeze-dried 
Pseudomonas  denilrifleans  (ATCC  13867)  were  revived  in  flasks  containing  125  mL  of 
nutrient  broth  (SIGMA.)  These  were  shaken  for  two  days,  after  which  plates  containing 
Tryptic  Soy  Agar  (SIGMA)  were  streaked  with  revived  cultures.  These  were  grown  in  an 
incubator  for  three  days  at  35°C  and  then  stored  at  4°C  in  a refrigerator. 

Growth  experiments  were  conducted  to  observe  behavior  upon  the  transition  from 
anoxic  to  aerobic,  as  well  as  aerobic  to  anoxic  growth  conditions.  Experiments  consisted 
of  three  phases.  A nutrient  solution  (Table  4-2)  containing  limited  carbon  substrate 
(0. 1 3 1 g/L  of  glutamic  acid)  was  used  for  the  initial  anoxic  phase.  The  solution 
contained  enough  carbon  for  bacteria  to  grow  to  an  approximate  absorbance  of  0. 1 , 
determined  by  calculation  and  confirmed  by  preliminary'  experiments.  This  phase  was 
carried  out  in  a reactor  (model  BioFlo  2000.  New  Brunswick  Scientific)  to  which  excess 
nitrate  (400  mg/L)  was  added.  Bacteria  were  added  to  the  reactor  directly  from  the  agar 
plates  and  grown  anoxically  at  30°C.  Nitrogen  was  bubbled  through  the  reactor  to  strip 
dissolved  oxygen. 

After  the  carbon  substrate  was  exhausted  in  this  anoxic  phase.  500  mL  of  the 
biomass  suspension  was  transferred  to  an  aerated  bioreactor  (model  F-1000.  New 
Brunswick  Scientific).  The  reactor  contained  nutrient  solution  with  excess  carbon  (5  g/L 
glutamic  acid)  to  allow  for  aerobic  growth.  The  rale  of  aeration  was  sufficient  so  that 
dissolved  oxygen  concentration  approached  air  saturation  in  the  reactor.  The  length  of 
this  aerobic  phase  was  varied. 
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The  final  anoxic  phase  began  after  the  end  of  the  aerobic  phase.  The  contents  of 
the  reactor  were  diluted  by  adding  nutrient  solution  containing  excess  carbon  substrate  (5 
g/L  glutamic  acid)  and  nitrate  (400  mg/L  as  N).  Nitrogen  gas  was  bubbled  through  the 
reactor  to  strip  dissolved  oxygen. 

An  automatic  sampling  system  was  utilized.  Individual  pumps  withdrew  culture 
media  from  the  reactor  and  passed  it  through  a flow  cell  inside  a spectrophotometer 
(Thermo-Spcctronic  Genesys  10UV).  The  spectrophotometer  was  programmed  to  take 
absorbance  measurements  at  550  nm  every  1 5 min.  throughout  all  three  phases. 


of  CuS04,  FeClj, 


Results  and  Discussion 

The  new  simpler  model  was  first  applied  to  the  results  of  experiments  performed 
by  Liu  et  al.  (1998b).  who  studied  the  effect  of  length  of  aerobic  phase  on  the  ensuing 

fit  was  obtained.  To  generate  these  fits,  the  same  parameters  were  used  for  both  runs 
(Table  4-3).  For  the  model  parameters  in  common  with  the  Liu  et  al.  (1998b)  model, 
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Fils  lo  two  experimental  runs  carried  out  in  the  present  research  are  shown  in 
Figure  4-2  and  4-3-  The  runs  differ  in  tile  length  of  the  intermediate  aerobic  phase  and 
concomitantly  the  lengths  of  the  diauxic  lags.  Parameters  and  unmeasured  initial 
conditions  used  for  the  fits  arc  given  in  Tabic  4-4.  A sample  Excel  spreadsheet  used  for 
the  fits  is  included  in  Appendix  A.  Figures  4-2  and  4-3  show  that  the  new  model  is 
successful  in  portraying  carbon-limited  anoxic  growth,  aerobic  growth  at  nonlimiting 
carbon  concentration,  and  anoxic  growth  at  nonlimiting  carbon  concentration.  The  model 
successfully  portrays  the  absence  of  lag  when  switching  from  anoxic  to  aerobic 
conditions.  It  also  successfully  portrays  the  length  of  the  diauxic  log  between  the  aerobic 
phase  and  the  final  anoxic  phase.  Furthermore,  the  model  is  successful  in  showing  the 
effect  of  aerobic  phase  length  on  the  subsequent  diauxic  lag  length.  The  experimental 
runs  shown  in  Figures  4-2  and  4-3  took  place  on  different  days.  Hence,  the  difference  in 
initial  enzyme  concentrations  reflected  in  the  table  would  be  expected.  All  of  the 
parameters  were  identical  with  the  exception  of  aml . Forcing  a common  value  on  that 
parameter  degraded  the  fits  somewhat. 
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Table  4-3:  Parameters  and  unmeasured  initial  conditions  for  fit  of  Liu  ct  al.  (1998b)  data 
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[ able  4-4:  Param 

eters  and  unntcasu 

red  initial  conditions  for  three-phase  experiments 

Parameters 

Figure  2 values 

Figure  3 values 

Po 

0.73 

0.73 

h" 

Pm 

0.336 

0.336 

h1 

Koh 

0.065 

0.065 

gO;  m ‘ 

Km 

0.77 

0.77 

gNOj-N  nr 

Ks, 

1 

1 

g COD  glutamic  acid  m ‘ 

K„m 

303.9 

303.9 

g COD  glutamic  acid  m " 

Kj 

0.5 

0.5 

gOnn5 

Yco 

0.9 

0.9 

g cell  COD/g  COD  glutamic  acid 

Yc.m 

0.9 

0.9 

g cell  COD/g  COD  glutamic  acid 

da, 

3.00  - 10"” 

3.00  • I0"" 

!>' 

am 

3.93  • 10-" 

3.93  - 10"“ 

h' 

am 

1.84  •10J)I 

1.84  • 10"" 

T1 

“mi 

5.38  • 10-01 

5.38  • 10"" 

h1 

" o 

0.5 

0.5 

”m 

3 

3 

b 

0.002 

0.002 

h-1 

Pn 

0 

0 

!>' 

Pm 

0.0002 

0.0002 

h1 

Unmeasured 

conditions 

Figure  2 values 

Figure  3 values 

Units 

Initial  e„ 

9.515  ■ I0J 

0 

g enzyme  COD/g  cell  COD 

Initial  eM 

1.375 

0.2367 

g enzyme  COD/g  cell  COD 

CHAPTERS 

EFFECT  OF  ENZYMES  AND  CARBON  SUBSTRATE  ON  THE  D1AUXIC  LAG 
It  has  been  suggested  that  the  pcriplasmic  nitrate  reductase  present  in  certain  bacteria 
plays  a role  in  their  transition  from  conditions  where  oxygen  is  present  to  conditions 
where  oxygen  is  absent.  To  that  effect,  the  diauxic  lags  of  denitrifying  bacteria 
containing  the  pcriplasmic  nitrate  reductase  ( Paracoccus  panloirophus.  formerly 
Thlospliaera  panlolropha)  were  compared  to  those  of  denitrifying  bacteria  lacking  the 
pcriplasmic  nitrate  reductase  ( Pseudomonas  dcnilrificans)  under  the  same  experimental 
conditions.  Furthermore,  since  one  of  the  primary  functions  of  the  pcriplasmic  nitrate 
reductase  is  redox  balancing,  a secondary  objective  was  to  investigate  the  influence  of  the 
oxidative  slate  of  the  carbon  substrate  on  the  transition  period  of  Nap-containing  bacteria. 

F.xpcrimcntal  Methods 

All  experiments  were  carried  out  under  aseptic  conditions.  Freeze-dried 
Paracoccus  panloirophus  (ATCC  355 1 2)  and  Pseudomonas  dcnilrificans  (ATCC  1 3867) 
were  separately  revived  in  flasks  containing  125  mL  of  nutrient  broth  (SIGMA.)  These 
were  placed  in  an  incubator  shaker  at  35°C  for  one  day  ( Pseudomonas  dcnilrificans 
flasks)  or  two  days  ( Paracoccus  panloirophus  flasks),  after  which  plates  containing 
Tryptic  Soy  Agar  (SIGMA)  were  streaked  with  revived  cultures.  These  were  placed  in 
an  incubator  for  three  days  at  35°C  and  then  stored  at  4°C  in  a refrigerator  for  no  more 


than  two  weeks. 
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Growth  experiments 

A premium:  phase  was  carried  out  before  each  experiment.  A culture  of 
Paracoccw  panlotrophus  or  Pseudomonas  denitrificam  was  prawn  aerobically  at  35°C 
in  nutrient  solutions  recommended  for  each  particular  culture:  ATCC  for  Paracoccus 
panlotrophus  (Table  5-1)  and  Komaros  el  al.  (1996)  for  Pseudomonas  denitrificam 
(Table  5-2).  Experiments  consisted  of  three  phases:  a startup  aerobic  phase,  a continuous 
aerobic  phase,  and  a batch  anoxic  phase.  All  three  phases  were  carried  out  in  a fermentor 
(BioFlo  2000.  New  Brunswick  Scientific). 

Nutrient  solution  containing  one  of  three  different  carbon  substrates  (malate. 
acetate,  or  caproate)  and  no  nitrate  was  used  for  the  startup.  A specific  amount  of 
biomass  from  the  preculture  was  added.  The  reactor  was  aerated  and  maintained  at  35°C. 
After  a short  period,  a feed  solution  containing  the  appropriate  carbon  substrate  at  a 
concentration  of  1000  mg/L  as  COD  (all  other  nutrients  in  concentrations  shown  in 
Tables  5-1  and  5-2)  was  pumped  into  the  reactor  to  yield  a dilution  rate  of  0.03h  '.  This 
dilution  rale  was  used  in  all  growth  experiments,  ensuring  diat  the  net  specific  growth 
rate  was  the  same  at  the  end  of  each  continuous  flow  period.  In  each  experiment,  the 
culture  was  grown  to  a steady  state  as  indicated  by  an  absorbance  (>.=550  nm.  I cm  path 
length)  variation  within  0,02%.  After  steady  state  was  achieved,  the  feed  was  stopped,  a 
portion  (typically  700-800  mL)  of  the  culture  was  rapidly  removed,  and  nutrient  solution 
containing  nitrate  (at  the  growth  temperature  of  35°C)  was  added  to  the  reactor,  giving  an 
in  situ  nitrate-nitrogen  concentration  of  at  least  200  mg/L.  Air  was  substituted  by 
nitrogen  gas  to  strip  out  dissolved  oxygen.  Growth  was  monitored  during  this  final  phase 
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by  taking  absorbance  measurements  al  550  nm  every  15  minutes.  Lags  and  specific 
growth  rates  were  determined  as  described  by  Lisbon  cl  a!.  (2002). 

Activities  of  Nap  and  Nar  were  determined  from  a biomass  sample  taken  near  the 
end  of  the  continuous  aerobic  phase.  Cells  were  harvested  by  centrifugation  (11,300 
RCF  for  10  minutes  at  4°C)  and  rinsed  with  20  mM  Tris  buffer  (pH  7).  They  were 
resuspended  in  2 mL  of  the  buffer.  The  assay  method  used  was  a modification  of  that 
given  by  Jones  ct  al.  (1977)  that  measures  the  activity  of  nitrate  reductase  in  intact  cells 
based  on  the  rate  at  which  reduced  benzyl  viologen  (BV)  decolorizes,  which  is 
proportional  to  the  rate  of  electron  transfer  to  nitrate.  Inside  an  anaerobic  chamber  (Coy 
Laboratories  Type  -A’),  several  3mm  glass  beads  were  added  (for  mixing)  to  a I cm 
optical  path  borosilicate  cuvette.  The  cuvette  was  filled  halfway  with  a solution  of 
benzyl  viologen  and  Tris  buffer  (0.3  mM  BV.  20  mM  Tris  HCI).  200  pL  of  resuspended 
biomass  were  added  (giving  a final  bacteria  concentration  of  approximately  55  mg/L  dry 
weight).  Enough  dilhionitc  (usually  25  pL  of  a 25  mM  solution)  was  added  to  reduce  the 
benzyl  viologen  to  give  an  absorbance  of  approximately  1 .8  at  550  nm.  More  benzyl 
viologenfTris  solution  was  added  to  completely  fill  the  headspace  of  the  cuvette.  The 
cuvette  was  then  stoppered  with  a Wheaton  seal  and  placed  in  a spectrophotometer 
(Thermo-Spectronic  Gcnesys  10UV)  programmed  to  take  automatic  measurements  every 
15  seconds.  After  monitoring  the  absorbance  for  3 minutes  to  give  a baseline,  the 
reaction  was  initiated  by  injecting  nitrate  through  the  rubber  seal  to  give  a final 
concentration  of  6 mM.  The  cuvette  was  inverted  twice  for  mixing  and  returned  to  the 
spectrophotometer,  where  the  absorbance  was  monitored  for  7 minutes.  The  seal  was 


then  broken,  air  was  allowed  into  the  cuvette  to  complete  the  oxidation  (decolorizaUon) 
of  benzyl  viologen,  and  the  remaining  absorbance  (absorbance  of  biomass  alone)  was 
measured.  The  slope  of  the  absorbance  vs.  time  plot  was  converted  to  units  of  benzyl 
viologen  oxidized  per  time,  which  would  be  equivalent  to  the  enzyme's  ability  to  catalyze 
the  use  reduction  of  nitrate.  The  unit  conversion  was  done  as  shown  in  Equation  5-1. 
Assuming  we  are  within  the  linear  range,  we  use 


(5-1) 


where  m is  the  slope  of  decolorizalion.  Cv  is  the  benzyl  viologen  concentration,  max  is 
the  maximum  absorbance  (which  includes  biomass  and  reduced  benzyl  viologen 
absorbance),  and  base  is  the  absorbance  of  biomass  only.  In  order  to  obtain  the  specific 
activity'  (concentration  change  per  biomass.  Equation  5-1  was  then  divided  by  biomass 
concentration  (in  units  of  mg/L  biomass). 


In  order  to  distinguish  between  activity  of  Nap  and  activity  of  Nar,  two  assays 
were  performed  simultaneously,  one  containing  I OOpM  azide  and  one  without  azide. 
Table  5-1 : Composition  of  nutrient  solution  for  hiracoccus  pamoirophus 
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Because  biomass  absorbance  is  directly  proportional  to  biomass  concentration,  an 
analogous  expression  applies  for  biomass  absorbance. 

A = -'•>  (5.5, 

\nA=(inA0-fim,la)+ftmm,i  (5-6) 

Here  A = biomass  absorbance,  in  = time  at  the  beginning  of  growth  in  the  aerobic  phase. 
An  = absorbance  at  the  beginning  of  the  aerobic  phase,  yw  = maximum  aerobic  specific 
growth  rate  (formerly  //„ ). 

Under  anoxic  conditions,  if  nitrate  is  in  excess,  as  is  the  case  at  the  beginning  of 
the  anoxic  phase,  an  analogous  equation  to  5-6  applies: 

lnA=(lnAl-i/„J.t,)  + Mm^,.l  (5-7) 

where  A = biomass  absorbance.  I;.  = time  at  end  of  diauxic  lag.  At  = average  absorbance 
during  the  lag.  maximum  anoxic  specific  growth  rate. 

In  order  to  calculate  the  maximum  specific  growth  rate  in  the  anoxic  phase.  In 
(absorbance)  vs.  lime  was  plotted.  Semi-logarithmic  plots  are  included  in  Appendix  B. 
A straight  line  is  fitted  and  the  slope,  according  to  Equation  5-7,  yields  the  maximum 
specific  growth  rate  for  anoxic  conditions.  The  end  of  the  lag  phase.  //.,  is  determined 
from  the  intersection  of  a line  of  zero  slope  drawn  at  the  overage  absorbance  during  the 
lag,  Al , and  the  straight  line  for  exponential  growth  at  the  beginning  of  the  anoxic  phase. 
The  diauxic  lag,  lo,  is  obtained  from 

'D-lL-ls  (5-8) 

where  is  is  the  time  at  which  air  was  turned  off. 
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In  some  experiments,  the  semi-logarithmic  plot  showed  that  during  the  second 
anoxic  phase  there  were  two  phases  of  exponential  growth,  each  with  different  growth 
rale.  In  these  cases,  two  straight  lines  were  fit  and  the  intersection  of  the  two  lines  was 
the  end  of  the  "lag". 

Results  and  Discussion 

Figure  5-1  presents  the  results  of  Pseudomonas  denilri/icans  growing  on  acetate. 
Tile  beginning  of  the  graph  shows  the  transition  from  the  startup  phase  to  the  continuous 
phase.  The  absorbance  then  becomes  approximately  constant  as  the  culture  approaches 
steady  state.  It  can  be  seen  that,  after  switching  from  the  aerobic  continuous  culture  to 
anoxic  batch  conditions,  the  culture  experienced  a diauxic  lag  of  16  hours  before 
resuming  exponential  growth.  Figure  5-2  shows  the  results  with  Paracoccus 
panioiroplius  growing  with  the  same  substrate.  As  before,  the  culture  was  switched  front 
aerobic  continuous  How  to  anoxic  batch  conditions  after  reaching  steady  slate.  An 
enlargement  of  the  anoxic  batch  phase  is  shown  within  Figure  5-2.  From  this  view  it  can 
be  seen  that,  after  switching  from  the  aerobic  continuous  culture  to  anoxic  batch 
conditions,  for  a period  of  approximately  1 hour  the  specific  growth  rate  was  relatively 
slow  (0.13  h'1),  after  which  it  increased  to  a higher  value  (0.21  h’1).  The  change  in 
growth  rate  is  easier  to  observe  in  the  semi-logarithmic  plot  in  Appendix  B (Figure  B-2). 
It  is  debatable  whether  or  not  this  should  be  considered  a diauxic  lag.  Traditionally  a 
diauxic  lag  has  been  defined  as  a period  of  little  or  no  growth  (Monod,  1942).  In  the 
present  case,  the  initial  growth  rate  is  62%  of  the  ultimate  growth  rate.  Whether  or  not 
this  is  to  be  called  a diauxic  lag,  the  period  over  which  the  rate  was  slower  is  short. 

Experiments  were  repeated  utilizing  a more  oxidized  carbon  substrate  (malate). 
Results  from  these  experiments  are  shown  in  Figures  5-3  and  5-4.  Once  again,  results 


diauxic  lag  (S  h)  whereas  i 


lag  (0.5 


In  order  to  verily  the  effect  of  carbon  substrate,  an  experiment  was  done  with 
Paracoccm  pantotrophus  using  caproatc  as  carbon  substrate.  Results  arc  shown  in 
Figure  5-5  and  the  semi-logarithmic  plot  is  included  in  Figure  B-5  in  Appendix  B.  After 

immediately  at  a constant  specific  growth  rate  (no  lag).  Pseudomonas  denitrificans  was 
unable  to  grow  using  caproale  as  carbon  substrate. 

shorter  lag.  followed  by  acetate  and  then  malate.  Clearly,  caproatc  had  no  lag.  and 
clearly  malate  had  some  lag.  albeit  a short  one.  As  mentioned  above,  it  is  debatable  as  to 

malate.  but  less  severe.  Table  5-3  shows  a summary  of  the  anoxic  specific  growth  rates 

in  Appendix  B. 

to  the  cell  membrane  in  intact  cells  (Bcdzyk  ct  al..  1999).  the  activities  measured  in  these 
assays  include  the  effect  of  Nar  and  Nap  combined.  Because  at  micromolar 
concentrations  azide  inhibits  Nar  and  has  no  effect  on  Nap.  this  inhibitor  can  be  used  to 

one  being  measured.  Because  oxygen  is  known  to  inhibit  Nar  and  these  samples  were 


45 

taken  after  a long  period  of  aerobic  growth,  no  Nar  activity  was  expected.  On 

assays  were  quite  high  compared  to  those  obtained  with  Pseudomonas  denilriflcans. 
which  were  practically  zero  (Table  5-4).  This  is  consistent  with  the  observed  longer  lags 

nature  of  the  carbon  substrate  used  for  growth. 

As  expected,  cultures  able  to  express  Nap  had  much  shorter  (if  any)  lags  than  those 
that  only  express  Nar  upon  switching  from  aerobic  to  anoxic  conditions,  thus  suggesting 
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Table  5-3:  Summary  of  experimental  results. 
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Table  5-5:  Parameters  and  unmeasured  initial  conditions. 


CHAPTER  6 
CONCLUSIONS 

A simplified  model  of  Ihc  diauxic  growth  of  bacteria  switching  between  electron 
acceptors  has  been  developed.  Based  on  the  application  of  this  model  to  new 
experimental  data  generated  as  pan  of  this  study,  as  well  as  data  from  the  literature,  the 
following  conclusions  may  be  drawn: 

• The  model  provides  excellent  (its  to  diauxic  growth  of  hcterouophic  bacteria 
switching  between  electron  acceptors. 


• The  model  is  capable  of  portraying  denitrification  dynamics  in  cyclic 
anoxic/aerobic  environments. 

• The  model  is  capable  of  simulating  carbon-limited  growth. 


Based  on  the  data  obtained  from  growth  experiments  using  denitrifying  bacteria 
containing  a penplusmic  nitrate  reductase  (Paracoccus pantolrophus ) and  denitrifying 
bacteria  containing  only  a membrane-bound  nitrate  reductase  {Pseudomonas 
denitrifwans),  the  following  conclusions  may  be  drawn: 

• The  pcriplasmic  nitrate  reductase  may  have  a role  in  speeding  the  transition  from 
growth  with  oxygen  as  terminal  electron  acceptor  to  growth  with  nitrate  as 
terminal  electron  acceptor. 


There  appears  to  be  a correlation  between  the  oxidation  state  of  the  carbon 
substrate  and  the  length  of  the  diauxic  lag. 


CHAPTER  7 
FUTURE  WORK 


Many  aspects  of  diauxic  behavior  when  switching  between  electron  acceptors  are 
still  a mystery.  Considerable  research  must  still  be  done  for  better  understanding  of  this 
phenomenon. 

One  key  question  yet  to  be  answered  is  how  the  carbon  substrate  concentration 
affects  the  length  of  the  diauxic  lag.  Experiments  described  in  Chapter  S where  cultures 
were  grown  continuously  under  aerobic  conditions  and  then  switched  to  batch  anoxic 
conditions  provided  low  substrate  concentrations  during  the  continuous  phase  and 
resulted  in  short  or  non-existent  diauxic  lags  for  Nap-containing  bacteria.  On  the  other 
hand,  experiments  not  shown  in  this  dissertation  where  Nap-containing  bacteria  were 
grown  under  aerobic,  batch,  excess  substrate  conditions  and  then  switched  to  batch 
anoxic  conditions  resulted  in  slightly  longer  lags  (as  long  as  3 hours).  Ellington  et  al. 


(2003)  studied  Nap  expression  in  continuous  cultures  of  Paracoccus  paMotrophus  in 
different  carbon  substrates  in  limited  concentrations  and  observed  that  levels  of 
expression  were  higher  than  those  in  batch  cultures.  Thus  it  seems  that  carbon  substrate 
concentration  inuy  have  an  effect  on  Nap  activity,  consequently  affecting  the  length  of 
the  lag.  One  way  to  examine  the  effect  of  carbon  substrate  concentration  is  to  perform 
continuous  experiments  with  different  dilution  rates. 

Our  observation  that  Nap  allows  for  bacteria  to  start  growing  anoxically 
immediately  after  aerobic  growth  leads  to  another  unanswered  question.  Microbiologists 
slate  that  because  nitrogen  reduction  using  Nap  is  not  directly  coupled  to  ATP 
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production,  it  docs  not  provide  a means  for  growth  and  is  only  used  for  redox  balancing. 
Therefore,  the  mechanism  for  growth  utilizing  the  periplasmic  enzyme  must  be 
examined.  The  fact  that  there  was  actual  growth  due  to  Nap  was  confirmed  by  batch 
experiments  where  Nap-containing  cultures  were  grown  anoxically  with  azide.  Because 
azide  inhibits  Nar.  the  growth  observed  during  these  experiments  was  attributed  to  Nap. 
Bedzyk  and  colleagues  ( 1 999)  also  observed  anoxic  growth  due  to  Nap.  In  fact, 
denitrification  by  Pseudomonas  sp.  strain  G-179  was  not  possible  without  Nap.  They 
could  not  define  the  mechanism  by  which  energy  is  gained  in  this  process.  A possible 
explanation  is  that  Nap  allows  for  bacteria  to  quickly  reduce  nitrate  to  nitrile,  and  the 
growth  observed  is  actually  growth  on  nitrite  instead  of  direct  growth  on  nitrate. 

However,  nitrile  reductase-deficient  mutants  of  Pseudomonas  sp.  strain  G- 1 79  can  grow 
on  nitrate  (Bedzyk  et  al.,  1999).  Regardless  of  these  findings,  the  hypothesis  should  be 
tested  with  Paracoccus paniotrophus  by  performing  growth  experiments  using  nitrite  as 
electron  acceptor  instead  of  nitrate.  If  long  lags  are  observed  afier  switching  from  using 
oxygen  to  using  nitrite  as  electron  acceptor,  this  mechanism  is  not  feasible. 

In  order  to  verify  that,  in  fact,  the  reason  why  Paracoccus  paniotrophus  has  no 
lag  is  because  it  contains  Nap.  it  is  important  to  perform  experiments  using  other  Nap- 
containing  bacteria.  Additionally,  it  would  be  ideal  to  try  to  inhibit  Nap  somehow  during 
growth  experiments  using  Paracoccus panioiroplius  or  to  cause  a mutation  in  this  species 
to  make  it  deficient  in  Nap  and  observe  if  this  leads  to  longer  lags. 

Because  the  application  of  this  study  is  in  wastewater  treatment  with  activated 
sludge,  another  suggestion  for  future  work  is  to  try  to  identify  what  percentage  of  bacteria 
in  activated  sludge  contains  Nap.  Our  research  group  is  currently  trying  to  isolate 
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bacteria  from  activated  sludge  for  this  purpose.  If  successful  in  isolating  cultures,  diauxic 
growth  experiments  should  be  performed  with  these  cultures. 


APPENDIX  A 

SAMPLE  EXCEL  SPREADSHEETS 
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=IF(F22»”(F22-H22JA2  "') 
=IF(F23<>“,(F23-H23)”2."-) 
=IF(F24<>~,(F24-H24)*2,'-) 
=IF(F25o_,(F25-H25)l'2.“) 
=IF(F26o-,(F26-H26)“2,“) 
=IF(F27<>-.(F27-H27)A2,“) 
=IF(F28<>”  (F28-H28)A2,“) 
=IF(F29<>“(F29-H29)A2,~) 
=IF(F30<>",(F30-H30)A2,-) 
=IF(F31  <>",<F31-H31  )A2,_) 
=IF(F32<>”(F32-H32)A2,-) 
=IF(F33»"(F33-H33)A2,-) 
=IF(F34o-',(F34-H34)A2,_) 
=IF(F35o"'.(F35-H35)A2,-) 
=IF(F36o"\(F36*H36)A2,"') 
=IF(F37o-.(F37-H37)A2,_) 
"IFfFSSo-.IFSB-HSS^,") 
=IF(F39<>"',(F39+I39)A2.-) 
=IF(F40<>-.(F40-H40)A2.") 
=IF(F4K»”.(F41-H41)A2,") 
=IF(F42<>“,(F42-H42)A2,-) 
=IF(F43<>-'.(F43-H43)A2,'"') 

=|F(F44«-,(F44-H44)A2~) 

=IF(F45<>“,(F45-H45)A2.“) 

=IF(F46<>“‘,(F46-H46)A2.'“) 

=IF(F47<>_,(F47-H47)A2,~) 

=IF(F48<>“,(F48-H48)A2,“| 

=IF(F49<>_,(F49-H49)A2.~I 

=IF(F50<>-,(F50-H60)A2,") 

=IF(F51<>~,(F51-H51)A2,“> 

=IF(F52<>”,(F52-H52)A2-') 

=IF(F53o-,(F53-H53r2.“) 

=IF(F64<>-,(F54-H54)A2,'-) 

=IF(F55<»~.(F55-HS5)A2,~) 
=IF(F86<>~,(F56-H56)A2,-) 
=IF(F57<>“(F57-H57)A2.-) 
=IF(FS8o~(F68-H58)A2,-) 
=IF(F59<>”.(F59-H59)A2.~) 
=IF(F60<>~,(F60-H60)A2,-) 
=IF(F61<>",(FB1-H61)A2,“) 
=IF(F62<>-'.(F62-H62)A2,-") 
=IF(F63<>"",(F63-H63)A2,-) 
=IF(F64-c>'_,(F64-H64)A2,“) 
= IFiF65-::-""(FS5-H65)A2,'-) 
=IF(F66o"(F68-H66)A2,"') 
=IF(F67o"",(F67-H67)A2,~) 
=IF(F68o"',(F68-H68)A2.-) 
=IF(F69o'"’,(F69-H69)A2.'"') 
=IF(F70o“.(F70-H70)A2.~) 
=IF(F71o“,(F71-H71)A2-') 
=IF(F72o-,(F72-H72)A2-') 


=H5O*Ot*T50 


=MAX<l22*OfU22.0) 
=MAX(l23*DfU23.0) 
=MAX(l24+OrU24,0) 
=MAX(l25*DrU25.0) 
=MAX(l26+Dt*U26,0) 
=MAX(I27+DI‘U27,0) 
=MAX(I28+DI'U28,0) 
=MAX(l29*Ot“U29,0) 
=MAX(I30*DIU30.0) 
=MAX(I31*D[-U31,0) 
=MAX(l32*DfU32,0) 
=MAX(l33*DfU33,0) 
=MAX(I34+DTU34,0) 
=MAX(l35*DtU35.0) 
=MAX(I36+DI*U36,0) 
=MAX(I37*DFU37,0) 
=MAX(I38*DFU38.0) 
=MAX(I39*DFU39,0) 
=MAX(l40*DfU40,0) 
=MAX(l41+Dt*U41,0) 
=MAX(l42*OfU42,0) 
=MAX(l43+Dt‘U43,0) 
=MAX(l44+Dt*U44,0) 
=MAX(I45*DFU45,0) 
=MAX(l46*DrU46,0) 
=MAX(I47+Dt'll47,0) 
=MAX(l48*DfU48.0) 
=MAX(I49A0!'U49,0) 
=MAX(l50»DfU50,0) 
=MAX(I51»DFU51,0) 
=MAX(l52*Dt’U52,0) 
=MAX(l53»DfU53.0) 
=MAX(l54+Dt'U64,0) 
=MAX(l55+DfU55.0) 
=MAX(l56*Dt'U56,0} 
=MAX(l57»DrU57.0) 
=MAX(l58+Dt*U58,0) 
=MAX(l59+Dt'U59,0) 
=MAX{l60+Dt*U60.0) 
=MAX(I61+DI*U61 .0) 

=MAX(I63*DFU63!o) 

=MAX(I64*DFU64,0) 

=MAX(I65+DFU65.0) 

=MAX(l66*Dt*U66,0) 

=MAX067+Dt'U67,O) 

=MAX(l68+Dt*U88,0) 

=MAX(l69*DfU69.0) 

=MAX(I70*DFU70,0) 

=MAX(I71*DI'U71,0) 


=MAX(J22+Dt*V22,0) 

=MAX(J23*Dt*V23,0) 

=MAX(J24+DrV24,0) 

=MAX(J25«DPV25,0) 

=MAX(J26+Dt*V26.0) 

=MAX(J28*DTV28]o) 

=MAX(J29+Dt*V29,0) 

=MAX(J30*DPV30,0) 

=MAX(J31*DfV31,0) 

=MAX{J32*DTV32,0) 

=MAX(J33*Dt*V33.0) 

=MAX(J34+Dt‘V34.0) 

=MAX(J35*DTV35,0) 

=MAX(J36*DrV36.0) 

=MAX(J37+Dt*V37,0) 

=MAX(J38*OTV38.0) 

=MAX(J39+Dt*V39,0) 

=MAX(J4<ROfV40,0) 

=MAX(J41*DrV41.0) 

=MAX(J43+Dt‘V43,0) 
=MAX(J44+Dt*V44.0) 
=MAX(J45+Dl‘V45,0) 
=MAX(J46»DrV46,0) 
=MAX(J47*Dt*V47.0) 
=MAX(J48*DrV48,0) 
=MAX(J49*CTV49,0) 
=MAX(J5(HDPV50,0) 
=MAX(  J51  +Dt*V51 .0) 

=MAX{J53*DrV53,0) 

=MAX(J54+01"V54,0) 

=MAX(J55*DrV55,0) 

=MAX(J5B*DTV56,0) 

=MAX(J57+DTV57.0) 

=MAX(J68*DrV58.0) 

=MAX(J59*Dt"V59,0) 

=MAX(J60*Dt*V60,0) 

=MAX{J61+Dt*V61,0) 

=MAX(J62*DrV62,0) 

=MAX(J63+Dt*V63.0) 

=MAX(J64+Dl*V64.0) 

=MAX(J65+Dt*V65.0> 

=MAX(J66+DPV66,0) 

=MAX(J67*D|-V67,0) 

=MAX(J68+DPV68.0) 

=MAX(J69*Di-V69.0) 

=MAX(J70*DfV70,0) 

=MAX(J71*Dn/71,0) 


=K25/H25 

=K26/H26 

=K28/H28 


=K3Q/H30 


=K38/H38 

=K39/H39 

=K40/H40 


=K46/H46 

=K48/H48 

=K49/H49 

=K51/H51 

=K52/H52 

=K53/H53 


=L26/H26 

=L30/H30 


=L43/H43 


=L49/H49 

=L50/H50 


-L52/H52 


=K59/H59  =159/H59  0 
=K60/H60  -L60/H60  0 
=K61/H61  =L61/H61  0 
=K62/H62  =L62/H62  0 

=K65/H65  =L65/H65  0 
=K66/H66  =L66/H66  0 
=K67/H67  =L67/H67  0 
=K68/H68  =L68/H68  0 
=K69/H69  =L69/H69  0 
=K70/H70  =L70/H70  0 
=K71/H71  -L71/H71  0 


=(30J*30_2-M22/eo_max)'022/(Koh*022ri22/(Ks_ox*l22)-H22 
H(aO_1+aO_2'M23yeo_max),023/(Koh+023)‘l23/(Ks_ox+l23)*H23 
=(aO_1+aO_2'M24/eo_max)‘024/(Koh+024)'l24/{Ks_ox+l24)‘H24 
=(aO_1*aO  2'M25/eo_max)’025/{Koh*025)‘l25/(Ks_ox+l25)’H25 
=(aO_1  *aO_2*M26/eo_max)*026/(Koh*026)*l26/(Ks_ox*l26)'H26 
s(aO_1+aO_2*M27/eo_max)'027/(Koh+027)*l27/(Ks_ox+l27)'H27 
=(aO_1+aO_2*M28/eo_max)'028/(Koh+028)'l28/(Ks_ox+l28)'H28 
=(aO  1+aO  2*M29/eo_max)'029/(Koh+029)'l29/(Ks  ox*I29)'H29 
=(aO_1  *aO_2*M30/eo_max)‘030/(Koh+030)'l30/(Ks_ox*l30)*H30 
=(aO_1+aO_2*M31/eo_max)'031/(Koh+031J‘l31/{Ks_ox+l31)*H31 
=(aO_1+aO_2‘M32/eo_max)‘032/(Koh+032>,l32/(Ks_ox+l32)‘H32 

=(aO_1*aO_2-M33/eo_max)'033/(Koh-033)'l33/(K5_ox-l33)'H33 
=(a0_1+a0_2'M34/eo_max)‘034/(Koh+O34)’l34/(Ks_0x+l34)’H34 
=(aO_1*aO_2*M35/eo_max)*035/(Koh+035)‘l35/(Ks_ox+l35)'H35 
=(aO_1+aO_2'M36feo_max)*036/(Koh+036)*l36/(Ks_ox+l36)*H36 
=(aO_1*aO_2-M37(co_max)-037/(Koh*037)-|37/(Ks_ox*l37)-H37 
=(aO_l*aO_2*M38/eo_max)*038y(Koh-‘038)'l38/(Ks_ox*l38|*H38 
=(aO_l*aO_2‘M39/eo_max)*039/(Koh*039}'!39/(Ks_ox*l39)'H39 
=(aCM*aO_2'M40/eo_max>'040/(Koll*040)'l40V(Ks_ox*MO)*H40 
a(aO_1+aO_2*M41/eo_max)'041/(Koh+041)'l41/(Ks_ox+l41),H41 
={aO_1+aO_2'M42/eo_max)*042/(Koh*042)M42/(Ks_ox+l42)'H42 
={aO_1*aO_2'M43/eo_max)*043/{Koh+043)*143/(Ks_ox*143)*H43 
=(aO_1*aO_2'M44/eo_max)-044/(Koh*044)-|44/(K5_ox*l44)'H44 
=(aO_1»aO_2'M45/eo_max)'O45/(K0h+O45n45/(K5_ox*l45)-H45 
=(aO_1+aO_2*M46/eo_max)'046/(Koh+046)*l46/(Ks_ox+l46)*H46 
=(aO_1+aO_2‘M47/eo_max)*047/(Koh+047)'l47/(Ks_ox*l47)‘M47 
=(aO_1+aO_2’M48/eo_max)'048/(Koh»048)*l48/(Ks~ox+l48)'H48 
=(aO_1*aO_2’M49/eo_max)*049/(Koh*049)'l49/(Ks_ox*l49)‘H49 
=(aO_1*aO_2-M50/eo_max)'050/(Koh-050)-|50/(Ks  OX*I50)'H50 
=(aO_l*aO_2*M51/eo_max)‘051/(Koh»051)'l51/(K5_oxx|5l)*H51 
=(aO_1+aO_2*M52/eo_max)’052/(Koh*052)‘IS2/(Ks_ox+l52)'H52 
=(aO_1*aO_2-U53/eo_max)-053/(Koh*053)-|53/(Ks_ox*l53)'H53 
=(aO_1»aO_2*M54/eo_max)*054/(Koh*054)*l54/(Ks_oxx|54)*H54 
=(aO_1+aO_2*M55/eo_max)*055/(Koh+055)'l55/(Ks_ox+l55)'H55 
=(aO_1+aO_2*M56/eo_max)’056/(Koh+056)"l56/(Ks_ox+l56)*H56 
=(aO_1+aO_2,M57/eo_max),057/(Koh+057J'l57/(Ks_ox+l57J'H57 
=(aO_1*aO  2*M58/eo  max)'058((Koh*058)-l58/(Ks_oxtl58)-H58 
=(aO_'*aO_2*M59/eo_max)‘059/(Koh+059)'l59/(Ks_ox+l59)*H59 
=(aO_1*aO_2*M60/eo_max)*060/(Koh*060)*l60/(Ks_ox»l60)‘H60 
=(aO_1*aO_2'M61/eo_max)'061/(Koh+061)'l61/(Ks_ox+l61)*H61 
=(aO_1  ♦aO_2‘M62/eo_max)'062/(Koh+062)‘l62/(Ks  ox*I62)'H62 

=(aO  1+aO  2-M63/eo_max)'063/(Koh*063)'l63/(Ks  0x»l63)'H63 
-(aO_1+aO_2*M64/eo_max|,064/(Koh+CD64}-l64/(Ks_ox+l64)’H64 
=(aO_1*aO_2,M65/eo_max|*065/(Koh+065)*l65/(Ks_ox+l65)’H65 
=(aO_1+aO_2*M66/eo_max)*066/(Koh+066)*l66/(Ks_ox+l66)'H66 
=(aO_1+aO_2*M67/eo_max)*067V(Koh+067)*l67/(Ks_ox+l67)*H67 
=(aO_l*aO_2'M68/eo_max)*068/(Koh+068)'l68/(Ks_ox+l68)*H68 
=(aO_1+aO_2'M69/eo_max)*069/(Koh+069)*169/(Ks_ox*l69)*H69 
=(aO_1+aO_2*M70/eo_max)'070/(Koh+070)*170/(Ks_ox*l70),H70 
=(aO_1*aO_2'M71/eo_max)'071/(Koh*071|'l71/(Ks_oxH71)'H71 
=(aO_1*aO_2‘M72/eo_max)*072y(Koh*072)*l72/(Ks~ox*l72)'H72 


=(aNO_HaNO_2‘N22/eno_max),J22/(Kno+J22p22/(Ks_an+l22)'H22‘Ktf(Ki+022) 

=(aNO_1+aN0_2'N23/ene_max)*J23/(Kno+J23)M23/(Ks_an+l23)'H23*Ki/(Kt*O23) 

=(aNO_1+aNO_2‘N24/eno_max)'J24/(Kno+J24)'l24/(Ks_an+l24)*H24'Ki/(Ki*024) 

={aNO_1*aNO_2*N25/eno_max}*J25/(Kno+J25)'l25/(Ks_an+l25)‘H25'Ki/(Ki+025) 

=(aNO_1+aNO_2‘N26/eno_max)M26/(Kfio+J26)'l26/(Ks_an+l26)*H26'Ki/(Ki+026) 

=(aNO_1+aNO_2’N27/eno_max)\)27/(Kno+J27)*l27/{Ks_an+l27|‘H27*Ki/{Ki+027) 

'(aNO_l*aNO_2-N28/eno_max)-J28/(Kno*J2e)-|28/(Ks_an*l28)4H28-Ki/(Ki*028) 

=(aNO_UaNO,2-N29/eno_max)'J29/(Kno»J29|'l29/(Ks_an»l29)-H29'Ki/(Ki.029) 

=(aNO_1+aNO_2,N30/eno_max)'J30/(Kno+J30)'l30/(Ks_an+l30)*H30,Ki/{Ki*030) 

={aNO_1+aNO_2*N31/eno_max(‘J31/(Kno+J31)'l31/{Ks_an+l31)-H31-Ki/(Ki+031) 

=(aNO_l*aN0_2*N32/eno_rnax)*J32/(Kno*J32)‘l32/(Ks_an*l32|‘H32‘Ki/(Ki*O32) 
={aNO_1*-aNO_2*N33/eno_max)‘J33/{Kno*J33),l33/{Ks_an*l33)‘H33'Ki/(Klx033) 
={aNO_1*aNO_2*N34/eno_max>‘J34/(Kno+J34)'l34/{Ks_an+l34)‘H34'Ki/(Ki+034) 
=(a^J*aNo“2-N3yeno_max)M35/(Kno^l35)'l35/(KsIan*l35)-H35,Ki/(Ki*035) 
=(aNO_1+aNO_2'N36/eno_max)*J36/(Kno+J36)-|36/(Ks_an+l36)*H36,KI/(KI+038) 
=(aNO_1+aNO_2’N37/eno_max>*J37/(Kno+J37)’l37/(Ks_an+l37)*H37‘Ki/(KiK>37) 
=(aNO_1*aNO_2*N38/eno_max}'J38/(Kno*J38)'l38/{Ks_an+l38)*H38*Ki/(Ki+038) 
={aNO_1*aNO_2"N39/eno_max}\l39/(Kno*J39)'l39/(Ks_an‘l39)‘H39‘Ki/(KL»039) 
=(aNO_1*aNO_2‘N40/eno_max)‘J40/(Kno*J40)‘l40/(Ks_an*l40)’H40‘Kl/(Kl*040) 
s{aNO_14aNO_2*N41/eno_max>*J41/{Kno+J41  )-|41/(Ks_an*l41  )'H4VKi/(Ki*041 ) 
=(aN0_!*aN0_2‘N42/eno_max)‘J42/(Kno*J42)‘l42/(Ks_an*l42)‘H42‘Ki/(Ki»042) 
=(aNO_1*aNO_2*N43/eno_max)*J43/{Kno+J43)‘l43/(Ks_an*l43)*H43*Ki/(Ki*043) 
=(aNO_1*aNO_2‘N44/eno_max)*J44/(Kno*J44)‘t44/(Ks_an*l44)*H44*Ki/(Ki*044) 
=(aNOJ-aNO_2-N45/enoImax)-J45/(Kno*J45)'l45/(KsIan*l45)-H45-Ki/(Ki*045) 
=(aNO_t*aNO_2*N46/eno_max)*J46/(Kno*J46)‘l46/(Ks_an*l46)'H46*Ki/(Ki*046) 
=(aNO_1*aNO_2-N47/eno_max)-J47/(Kno*J47)-|47((Ks_an*l47)-H47-Kl/(Ki*047) 
=(aNO_1*aNO_2*N48/eno_max)'J48/(Kno*J48)*l487{Ks_an*l48)‘H48'Kl/(Ki*048) 
={aNO_1*aNO_2‘N49/eno_max)’J49/(Kno»J49)*l49/(Ks_an*l49),H49‘Kl/{Ki*049) 
=(aNO_l*aNO_2-N50/eno_max)-J5(V(Kno*J50)'l50/(Ks_an*l50CH50-Ki/(Ki*050) 
s{aN0_1+aN0_2'N51/eno_max)M51/(Kno*J51)’l51/(Ks_an*l51)'H51‘Ki/(Ki+051) 
=(aN0_1*aN0_2‘N52/eno_max)‘J52/(Kno*J52)*l52/(Ks_an*l52)‘H52'Kl/(Ki*052) 
=(aNO_1*aNO_2‘N53/eno_max|’J53/(Kno*J53)'l53/(Ks_an*l53)'H53‘Ki/(Ki*053) 
=(aNO_1-*aNO_2*N54/eno_max)*J54/(Kno+J54)‘l54/(Ks_an+l54)‘H54,Ki/(Ki*054) 
=(aNO_1+aNO_2'N55/eno_max)*J55/(Kno*J55)'l55/(Ks_an+l55)*H55'Ki/(Ki+055) 
=(aNO_1+aNO_2'N56/eno_max)*J56/{Kno*J56)‘l56/{Ks_an+l56)*H56‘Ki/(Ki*056) 
=(aNO_naNO_2-N57/eno_max)\l57/(Kno*J57)-l57/(Ks_an*l57),H57‘Ki7(Kl*057) 
=(aN0_1*aN0_2'N58/eno_max)*J587(Kno*J58)'l587(Ks_an*l58f,H58*Ki/(Ki*O58) 
=(aNO_1*aNO_2-N59/eno_maxCJ597(Kno*J59),l597(Ks_an*l69rH59-Ki/(Ki*059) 
={aNO_1*aNO_2'N60/eno_max)‘J60/(Kno*J60)'l60/(Ks_an*l60)'H60‘Ki/(Ki*060) 
=(aNO_1+aNO_2‘N61/eno_max)"J61/(Kno*J61),l61/{Ks_an+l61)*H6rKi/(KI*061) 
=(aNO_l»aNO_2'N62/eno_max)'J62/(Kno*J62)*l62/(Ks_an*l62)"H62’Ki/(Ki*062) 
=(aNO_1  *aNO_2‘N63/eno_max)-J63/(Kno*J63)*l63/(Ks_an*l63)'H63‘K[/(Ki*063) 
=(aNO_1*aNO_2,N64/eno_max)‘J64/(Kno*J64)’l64/(Ks_anH64)'H64-KI/(Ki»064) 
=(aNCM  *aNO_2'N65/eno_max)*J65/(Kno*J65)‘l65/(Ks~an*l65)'H65‘Ki/(Kl*065) 
=(aNO_1*aNO_2,N66/eno_max)"J66/(Kno+J66)*l66/(Ks_an+l66)*H66,Ki/(Ki+066) 
=(aNO_1+aNO_2'N67/eno_max)*J67/(Kno+J67)*l67/(Ks_an+167)*H67*Ki/(KI+067) 
=(aNO_1+aNO_2'N68/eno_max)\l68/(Kno+J68)*l68/(Ks_an+l68)*H68*Ki/(Ki+068) 

=(aNO_1*aNO_2-N69/eno_max)-J69/(Kno*J69)-|69/(Ks_an*l69|-H69-Ki/(Ki*069) 

=(aNO_1»aNO_2*N70yeno_max)’J70/(Kno*J70)*l70/(Ks_an*t70)‘H70*Ki/(Kl*070) 

={aNO_l*aNO_2'N71/eno_max)‘J71/(KnoxJ71)*l71/(Ks_an«l71)*H7TKi/jKi*071) 

=(aNO_1*aNO_2‘N72/eno_max}’J72/(Kno*J72)*l72/(Ks_an»l72)'H72‘Ki/(KI*072) 


=mumax_ox‘(M22/eo_max)An_ox*022/(Koh+022)'l22/(Ks_<w 

=mumax_ox‘(M23/eo_max)An_ox*023/(Koh+023)'l23/(Ks_ox 

=mumax_ox-(M24/eo_max)”n_ox-02a/(Koh*024)-|24/(Ks_o. 

=mumax_ox‘{M26/eo_max)"n_ox*025/(Koh*025)-l25/(Ks_o> 

-mumax_ox‘(M26/eo_max)'n_ox'026/(Koli*026)*l26/(Ks_o« 

=rTHimax_ox*(M27/eo_max)'n_ox’027/(Koh*027)*l27/(Ks_o> 

-mumax_ox*(M28/eo_max)"n_ox‘028/(Koh*028)'l28/(Ks_o> 

=mumax_ox'(M29/eo_max)"n_ox‘029/(Koti*029)-|29/(Ks_ox 

=mumax_ox*(M30/eo_max)"n_ox'030/(Koh*030)'l30y(Ks_o» 

=mumaxIox'(M31/eo_max)”n~ox-031/(Koti*031)-l31/(Ks_o> 

=mumax_ox'(M32/eo_maxlAn_ox*032/(Koh+032)*l32/(Ks_o> 

=mumax_ox'(M33(eo_max)'n_ox-033/(Koti*033)-|33/(Ks_o> 
=mumax_ox'(M34feo_max)An_ox'034/(Koh+034)*l34/(Ks_o* 
=mumax_ox'(M35/eo_max)"n_ox-035/(Koh*035)'l35/(Ks_o> 
=mumax_ox-(M36/eo_max)"n_ox'036/(Koti*036)'l36/(Ks_o> 
=mumax_ox‘(M37/eo_max)"n_ox‘037/(Koh*037)*l37/(Ks_o> 
=mumax_ox'(M38/eo_max)*n_ox-038/(Koh*038)-|38/(Ks_o. 
=mumax_ox-(M39/eo_max)“n_ox-039/(Koh*039)'l39/(Ks_o> 
=mumax_ox‘(M40/eo_max>An_ox*040/(Koh*040)'l40/(Ks_o> 
-mumax_ox*(M41/eo_max)"n_ox*041/(Koh*CW1|*141/(K$_o» 
=mumax_ox'(M42/eo  max)»n  ox'042/(Koh*042)-|42/(Ks_o» 
=mumax_ox-(M43/eo_max)*n_ox,043/(Koh*043)-l43/(Ks_o> 
-mumax_ox*(M44/eo_max)An_ox‘044/(Koti*044)M44/(Ks_ox 
=mumax_ox,(M45/eo_max>An_ox*045/(Koh+045)*l45/(Ks_o> 
=mumax_ox'(M46/eo_max)An_ox*046/(Koh+046)'l46/(Ks_o> 
=mumax_ox*(M47/eo  max)An_ox'C347/(Kon*047)-l47/(Ks_o. 
=mumax_ox*(M48/eo_max)An_ox'048/(Koh+04B)*M8/(Ks_o> 
=mumax_ox’(M49/eo_max)An_ox*049/(Koh+049)*l49/(Ks_ox 
=mumax_ox’(M50/eo_max)An_ox*050/(Koh*050)’l50/(K$_oi 
=mumax_ox*(M51/eo  max)An_ox‘05 1 /(Koh*05 1 )‘l  5 1 /(Ks_o> 


)An_ox‘053/(Koh*053)'l53/(Ks_o« 
ix_ox*{M54/eo_max)An_ox'054/(Koh+054)*l54/(Ks_ox 
ix_ox*(M55/eo_max)An_ox'055/(Koh+055)'155/(Ks_o> 
ix_ox*{M56/eo_max)An_ox'056/(Koh+056)*l56/(Ks_o> 
lx_ox'{M57/eo_max)An_ox*057/(Koh+057),l57/(Ks_o> 
ix_ox'(M58/eo_max)An”ox,058/(Koh+058)*l58/(Ks_o> 
iax)An_ox*059/(Koh*059)*l59/(Ks_o< 


lx_ox'(M63/eo_max)An_ox‘063/(Koh+063)"l63/(Ks_oxAl63)*H63 
lx_ox*{M64/eo_max)An_ox‘064/(Koh+064}*l64/(Ks_oxA164)'H64 
lx_ox-(M65/eo_max)An_ox'066/(Koh*065)'l66/(Ks_ox*l65)'H65 
ix_ox'(M66/eo_maxiAn_ox*066/(Koh+066)'l66/(Ks_ox+l66)'H66 
ix_ox'(MB7/eo_max)An_ox'067/(Koh+067)'l67/(Ks_ox»t67)*H67 
<S_0X*I68)*H68 
(s_ox*I69)‘H69 
An_ox’07Q/(Koh+070)*l70/(Ks_ox+l70)*H70 
ix_ox*(M71/eo_max)An_ox'071/(Koh+071)‘l71/(Ks_ox+l71)'H71 
lx_ox’(M72/eo”max)An~0X‘O72/(Koh+O72)*l72/(Ks~ox*l72)'H72 
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_an‘(N22/eno_max)An_anox‘J22/(Kno*J22)'l22/(Ks_arv*l22)*H22'Ki/(Ki*022) 
_an*(N23/eno_max)An_anoxM23/{Kno+J23)'l23/(Ks_atvH23)‘H23'Ki/(Ki+023) 
_an*(N24/eno_max)An_anox\J24/{Kno+J24)'l24/(Ks_an+l24),H24'Ki/(Ki+024) 
_an‘(N25/ei>0_max)An_anox\J25/(Kno*J25)‘l25/(Ks_an+l25)'H25‘Ki/(Ki+O25) 
_an’{N26/eno_max)Rn_anox*J26/(Kno+J26)-|26/(Ks_an+l26)*H26*Ki/(Ki+026) 
_an-(N27/eno_max)Rn_anox'J27/(Kno*J27),l27/(Ks_an*l27)'H27-Ki/(Ki*027) 
_an*(N28/eno  max)An  anox’J28/(Kno*J28)-|28/(Ks_an*!28)'H284Ki/(Ki*028) 
_an*(N29/eno_maxJAn_anox‘J29/(Kno*J29)*l29/(Ks_an+l29)’H29*Ki/(KH029) 
_an*(N30/eno_max)An_anox*J3Q/(Kno*J30)*l30/(Ks_an*l30)‘H30*Kl/(Ki*030| 
_an-(N31/eno_max)"n_anox-J31((Kno*J31)-l31/(Ks_an*l31)-H31-Ki/(Ki*031) 
_an*(N32/eno  max)An  anox-J32/(Kiro»J32)'l32/(Ks  an*l32)'H32-Ki/(Ki*032) 
_an'(N33/eno_max)An_anox‘J33/(Kno*J33)'l33/(Ks_an+I33)‘H33'Ki/(KH033) 
_an*(N34/eno_max)"n_anox'J34/(Kno*J34)‘134/(Ks_an*l34)'H34'Kl/(Ki*034) 
_an'(N35/eno_max)"n_anox‘J35/(Kno*J35)'l35/(K5_an*l35)'H354Ki/(Ki*035) 
_an'(N36/eno  max)"n  anox'J3»(Kno*J36)-|36/(Ks  an*l36)‘H36'Ki/(Ki-036) 
=mumax_an*(N37/eno_max)An_anox-J37/(Kno*J37)*l37/(Ks_an+l37),H37'KI/(KH037) 
_an'(N38/eno_max)An_anox'J38/(Kno4J38)‘l38/(Ks_an+l38)'H38'Ki/(Ki-*038) 
_an*(N39/eno_max)An_anox‘J39/(Kno+J39)‘l39/(Ks_an+l39)*H39*Ki/(Ki*039) 
• _an’(N40/eno_max)An_anox'J40/(Kno*J40)*l40/(Ks_an+l40)'H40*Ki/(Ki*040) 
_an'(N41/eno_max)An_anox,J41/(Kno+J41)4M1/(Ks_an+l41),H41*Ki/(Ki+041) 
_an*(N42/eno_max)An_anox'J42/(Kno+J42)*l42/(Ka_an+l42)'H42'KI/(Ki+042) 

_an-(N43/eno_max)An_anox*J43/(Kno*J43)'l43/(Ks_an*l43)*H43'Ki/(KI*043) 
_an'(N44/eno_max)An_anox,J44/(Kno+J44)4144/(Ks_an+l44)*H44'Ki/(Ki+044) 
'_an*(N45/eno_max)An_anox‘J45/(Kno+J45)*l45/(Ks_an+l45)‘H45'Ki/(Ki+045) 
_an'(N46/eno_max)An_anox‘J46/(Kno+J46)-l46/(Ks_an+146)'H464Ki/(Ki+046) 
_an’(N47/eno_max)An_anox*J47/{Kno*J47),l47/(Ks  an*l47)'H47'Ki/(KI*047) 
_an'(N48/eno_max)An_anox‘J48/(Kno+J4B)'l48/(Ks_an+l48)'H48‘Ki/(Ki+048) 
_an*(N49/eno_max)An_anox*J49/(Kno+J49)AM9/(Ks”an+l49)*H49,Ki/(Ki+049) 
_an4(N50/eno_max)An_anox*J60/(Kno+J50)4l50/(Ks_an+l50)'H50'Ki/(Ki+050) 
_an*(N51/eno_max)An_anox\l51/(Kno+J51)*l51/(Ks_an+l51)*H51‘Ki/(Kj+051) 
_an*(N52/eno_max}An_anox'J52/(Kno+J52)*l52/(Ks_an+l52)'H52'Ki/(Ki+052} 
_an*(N53/eno_max)An_anox*J53/(KnoAJ53)‘l53/{Ks_an+l53)AH53‘Ki/(Kt+053) 
_an-(N54/eno_max)An_anox-J54/(Kno*J54)'l54/(Ks_an*l54)'H54-Ki/(Ki*054) 
_an*(N55/eno_max>An_anoxM55/(Kno*J55)’l55/{Ks_an+l55)‘H55'Ki/(Ki+055) 
_an,(N58/eno_max|An_anox*J56/(KnoAJ56)'l56/(Ks_an+l56)*H56AKi/(Ki*056) 
_an'(N57/eno_max)An_anox\l57/(Kno*J57)*l57/(Ks_an+l57)*H57*Ki/(Ki+057) 
_an*(N58/eno_max)An_anox*J58/(Kno+J58)‘l58/(Ks_an+l58)‘H58'Ki/{Ki*058) 

_an*(N59/eno_max)An_anox*J59/(Kno*J59)*l59/(Ks_an*t59)’H59*Kt/(Ki»059) 
_an-(N60/eno_max)An_anox'J60/(Kno«J60)'l60((Ks_an*l60CH604KI/(KI*060) 
_an4(N61/eno_max)An_anox'J61/(Kno*J61)-|61/(Ks_an*l61)-H6rKL'(Ki*061] 
_an'(N62/eno_max)An_anox'J62/(Kno*J62)‘l62/(Ks_an*t62)*H62‘Ki/(Kt*062) 
_an*(N63/eno_max)An_anox'J63/(Kno*J63)‘!63/(Ks_an*l63)'H63'Ki/(Kl*063) 
_an'(N64/eno_max}An_anox'J64/(Kno+J64)‘l64/(Ks_an*l64)‘H64'Ki/(Ki+064) 
_an,(N65/eno_max)An_anox’J6S/(Kno+J65)*l65/(Ks_an+l65)*H65’Ki/(KI*065) 
_an4(N66/eno_max)An_anox\l66/(Kno+J66)*l66/(Ks_an*l66)'H66*Ki/(Ki+066) 
_an'(N67/eno_max)An_anox\l67/(Kno*J67)'l67/(Ks_an+l67)*H67'Ki/(Ki+067) 
_an‘(N68/eno_max)An_anox'J68/(Kno+J68>*l68'0<s_an+l68)'H68*Ki/(Ki+O68) 
_an'(N69/enQ_ntaxJAn_anox'J69/(Kno*J69)‘l69/(Ks_an-‘l69)*H69’Ki/(Ki*069) 
_an’(N70/eno_max)An_anox4J70/(Kno+J70)170/(Ks  an*l70)'H70'KI/(Ki«070) 
_an*(N71/eno_max)An_anox‘J71/(Kno+J71)*l71/(Ks_an+l71),H71*Ki/(Ki+071) 
_an'(N72/eno_max)An_anox'J72/(Kno+J72}'l72/(Ks_3n+l72)'H72*Ki/(Ki*072) 
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=-((l/Yc_an)/{2  86*Yc_an))*S22 
=-<(Wc_an)((2  86‘Yc_an))*S23 
=-<(1'Yc_an)/(2.86'Yc_an))-S24 
H(t/yc_an)/(2  86-Yc_an))-S25 

=-((1/Yc~an)/(2.86"Yc~an))'S27 
H(1/Yc_anW2.86-Yc_an))-S28 
HIl'Yc  an)/(2  86*Yc  an))‘S29 

=-f(1/Yc_an)/(2  86*Yc_an))*S30 
HC'Yc_an)/(2  86-Yc_an))-S31 
=-<(1/Yc_an)/(2  66*Yc_an))*S32 
=-|(1lYc_anW286-Yc_an))-S33 
=-«  Wc_an)/(2  86-Yo_an))‘S34 
H(1/Yc_an)((2.86-Yc_an))-S35 

H(1/Yc  an)/(2.86*Yc  an))'S36 
=-«1/Yc_anW2.86-Yc_an)rS37 
H(Wc_an)((2.86-Yc_an))-S38 
=-«1/yc_an)/(2.86-Yc_an))-S39 
H(1/Yc_an)f(2  86-Yc_an))-S40 
H(1/Yc_anV(2  86-Yc.an))-S41 

= ((17YcIan)/(286-YcIan))-S« 
H(1'Yc_an)/(2  86-Yc_an))‘S44 
=-<(1/Yo  an)/(2  86"Yc  an))‘S45 
=-((VYc_ar)/(2.86'Yc  an))‘S46 
=-((1'Yc_anW2.86-Yc_an))'S47 
H(1'Yc_an)fl2.86-Yc_an))-S48 
=-((1/Yc_an)/(2  86'Yc_an))‘S49 
HC'Yc  an)/(2  86‘Yc  an))'S50 
=-<(1/Yc_anl/(2  86*Yc_an))*S51 
=-((’"0  an)/(286*Yc  an))'S52 
=-«1A’c_an)/(2  86-Yc_an))‘S53 
=-((''Yc_an)/(2  86-Yc_an))-S54 
=*<(1^Yc_an)/(2.86*Yc_an))*S55 
H(Wc_an)/(2.86-Yc_an))-S56 
=-((1/Yc_an)/(2.86-Yc  an))-S57 
H(1'Yc_an)/(2.86-Yc_an))-S58 

=-«17Yc_an)/(2  86-Yc_an))-S59 

=-((1/Yc~an)/(2.86'Yc”an))'S61 
=-H1^Yc_an)/(2.86*Yc_an))*S62 
=-((1/Yc  an)/(2.86*Yc_an))‘S63 
= ((1/Yc_an)/(2  86-Yc_an))-S64 
=-«  Wc_an)/(2  86-Yc_an))‘S65 
=-«1'Yc_an)/(2  86-Yc_an))-S86 
=-((17Yc_an))(2  86-Yc_an))>S67 
H(1/Yo  an)/(2  86*Yc  an))'S68 
H(1'Yc_anW2.86-Yc.an))-S69 
H(Wc_an)«2.86-Yc_an)rS70 
H(Wc_anW2.86-Yc_an))*S71 
=-((1A'c_an)/(2.86-Yc_an))'S72 
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=IF(F25<>~,(F25-H25)A2."") 
=IF(F26<>-.(F26-H26)"2."') 
=IF(F27<>”,(F27-H27)»2.“) 
=IF(F28o~,(F28-H28)»2,") 
=IF(F29<>”,(F29-H29)A2  ™) 
=IF(F30<>~,(F30-H30)A2,'"') 
=IF(F31  o-.(F31-H31  y2."") 
=IF(F32*>-.(F32-H32)A2.”) 
=IF(F33<>_,(F33-H33)*2.*7 
=IF(F34»-(F34-H34)A2.'-') 

=IF(F35<»_,(F35-H35)'2.-) 

=IF(F38<>”.(F36-H36)*2,") 
=IF(F37<>-.(F37-H37)-2.~) 
=IF(F38<>”,(F38-H38r2,~) 
=IF(F39<>-.(F39-H39)»2.-) 
=1F(F40<>-.(F40-H40)*2,”) 
=IF(F41  <>“,(F41-H41)-2.“) 
=IF(F42<>“.(F42-H42)A2 ,”) 
=IF(F43<>~,(F43-H43)“2,-) 
=IF(F44o~(F44-H44)»2,”) 
=IF(F4S<>“.(F46-H45)»2,") 
=IF(F46<>'",(F46-H46)*2,-) 

=IF(F47<>-,(F47-H47)‘2.’-) 

=IF(F48<>“,(F48-H4B)“2,-) 

=IF(F49<>~(F49-H49)*2.-) 

=IF{F50<>“.(F50-H5Q)*2,-) 

=IF(F51«~,(F51-H51)*2,-) 

=IF(F52<>”,(F52-H52)»2,-) 

=IF(F53<>”",(F53-H53)*2,_) 

=IF(F54<>~(F54-H54)»2,-) 

=IF(F55<>"'.(F55-H55)*2,~) 

=IF(F56»-'.(F56-H56)*2,_) 

=IF(F57<>"'.(F57-H57)»2.“) 

=IF(FS8<>_.(F58-H68)”2,-) 

=IF(F59<>“,(F59-H59)A2,") 

=IF(F60o-,(F60-H60)»2.-) 

aIF(F61o"",(F61-H61)A2,”) 

=IF(F62<>"'.(F62-H62)A2,“) 

=IF(F63<>"',(F63-H63)A2,“) 

=IF(F64o-',(F64-H64)A2,“) 

=IF(F65<>~(F65-H65)A2,-) 

=IF(F66<>“.(F66-H68)A2,-) 

=IF(F67o-'.(F67-H67)A2.“) 

=IF(F68<>"".(F68-H68)A2.") 

=IF(F69<>_.(F69-He9)A2,”) 

=IF(F70-o'-.(F70-H70)A2,-> 

=IF(F7K>-,(F71-H71)A2.~) 

=IF(F72<>".(F72-H72)A2,“) 

=IF(F73o-",(F73-H73)A2,~) 

=IF(F74<>'",(F74-H74)A2,'") 

=IF(F75o",(F75-H75)A2,“) 


=H59+DtT59 

-H60*D1*T60 

=H62*DIT62 


=MAX(J25+DrV25.0) 
=MAX{J26*Dt‘V26,0) 
=MAX(J27*D1-V27.0) 
=MAX(J28+DTV28.0) 
=MAX(J29*DfV29,0) 
=MAX(J30+Dt*V30,0) 
=MAX(J31*DTV31.0) 
=MAX(J32+Dt*V32.0) 
=MAX(J33+Dt*V33,0) 
=MAX(J34*DfV34.0) 
-MAX(J35+Dl"V35,0) 
sMAX(J36+Dt*V36,0} 
=MAX(J37*DfV37,0) 
=MAX(J3B+Dl*V38,0) 
=MAX(J39*DTV39.0) 
=MAX(J40+Dt*V40.0) 
=MAX(J41+DrV41 ,0) 
=MAX(J42+Dt*V42,0) 
=MAX(J43*DfV43,0) 
=MAX(J44*D|-V44,0) 
=MAX(J45*DI-V45,0) 
=MAX(J46*DfV46,0) 
=MAX(J47+DI*V47,0) 
=MAX(J48+Dl*V48.0) 
=MAX(J49»OfV49,0) 
=MAX(J50K3TV50.0) 
=MAX(J51+DrV51,0J 
=MAX(J52+DI"V52,0) 
=MAX(J53»DrV53,0) 

=MAX(J54*DfV54,0) 

=MAX(J55+Dt*V55,0) 

=MAX(JS6*DrV56.0) 

=MAX(J67+OrV57,0) 

=MAX(J58*DTV58.0| 

=MAX(J59*Dl*V59.0) 

=MAX(J60*DrV60,0) 

=MAX{J61+Dt*V61.0) 

=MAX<J62*DfV62.0) 

=MAX(J63*DI-V63,0| 

=MAX(J65+DTV6S.O) 

=MAX(J66+DrV86.0) 

=MAX(J67+Dt*V67.0) 

=MAX(J68+DI*V68,0) 

=MAX(J69+Dt*V69.0) 

=MAX(J70*Dt-V70.0) 

=MAX{J71+Dt*V71,0) 

=MAX{J72+Dl*V72,0) 

=MAX(J73*DI'V73,0) 

=MAX(J74+Dt’V74,0) 


=K49+Dl*W49 

=K50*DfW50 

=K53*DfW53 

=K55*D1*W55 


sUssSIsiii 


=L26/H26 

=L30/H30 

=L33/H33 

=L38/H38 

=L40/H40 


=L48/H48 

=L58/H58 

=L60/H60 

=L62/H62 

=163/H63 

=L65/H65 


=L68/H6B 

=L69/H69 
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=(aO_I»aO_2*M25/eo_max)‘025/(Koh*025)‘l25/(Ks_OX*l25)‘H25 

=(aO_1+aO_2*M26/eo_max)'02G/(Koh*026)’l26/(Ks_OX*l26)'H26 

={aO_1+aO_2*M27/eo_max)'027/(Koh+027)*l27/(Ks_ox+l27)'H27 

={aO_1+aO_2'M28/eo_max)‘028/(Koh*028)*l28/(Ks_ox+l28)‘H28 

=(aO_1+aO_24M29/eo_max)4029/(Koh+029)4l29/(Ks_ox+l29)4H29 

s(aO_1+aO_2*M30/eo_max)‘03Q/(Koh*030)'l30/(Ks_ox*l30)*H30 

=(aO_1+aO_24M31/eo_max)4031/(Koh+031)'l31/{Ks_ox+l31)4H31 

=(aO_1+aO_2'M32/eo_max)'032/(Koh-*032)'l32/(Ks_ox+l32)'H32 

=(aO_1+aO_2*M33/eo_max)*033/(Koh+033)*l33/(Ks_ox+l33)*H33 

-(aO~1*aO~2‘M34/e0~max)"O34/(K0h*O3d)*l34/(Ks_ox*l34)*H34 

=(aO_1+aO_24M35/eo_max)4036/(Koh+035)4l35/(Ks_ox+l35)4H35 
=(aO_1+aO_2*M36/eo_max)‘036/(Koh*036)'l36/(Ks_ox+l36)'H36 
=(aO_1*aO_24M37/eo_max)‘037/(Koh*037)-l37/(Ks_OX*l37)'H37 
=(aO_1*aO_2'M38/eo_max)-038/(Koh*038)4l38/(Ks_ox*l38)'H38 
={aO_1+aO_2*M39/eo_max)4039/(Koh+039)4l39/(Ks_ox+l39)4H39 
=(aO  1*aO_2*M40/eo  max)4040/(Koh»040)4l40/(Ks_ox»l40)4H40 
=(aO_1+aO_24M41/e0_max)4O41/(K0h*O41)4l41/(Ks_0X*l41)4H41 
=(aO_1*aO_24M42/eojlu»)4CM2/(Koh*042)4l42/(Ksox*l42)4H42 
s(aO_1+aO_24M43/eo_max)4043/(Koh+043)4l43/(Ks_ox+l43)4H43 

=(aO_l*aO_2-M44/e0_max!'O44/(K0h*O44)'l44/(Ks_ox»l44)'H44 

=(aO_1*aO_2-M45/eo_max|'045/(Koll*045)'l45/(Ks_ox*l45)-H45 

=(aO_l*aO_2‘M46/eo_max)-O46/(K0h*O46)'l46/(Ks_ox*l46)-H46 

=(aO_1*aO_24M47/eo_max)'047/(Koh*047)'l47/{KS_OX*l47)'H47 

=(aO_1+aO_24M48/eo_max)"048/(Koh+048)4l48/(Ks_ox+l48)4H48 

=(aO_1*aO_24M49/eo_max)4049/(Koll*049)4l49/(Ks_ox»l49)4H49 

=(aO_1+aO_24M50/eo_max)4050/(Koh*050)4l50/{Ks_ox*l50)4H50 

=(aO_1+aO_24M51/eo_max)4051/(Koh+051)"l51/(Ks_ox+l51)4H51 

=(aO_1+aO_24M52/eo_max)4052/(Koh*052)4l52/(Ks_ox+l52)"H52 

=(aO_1+aO_24M53/eo_max)4053/(Koh+053)4l53/(Ks_ox+l53)4H53 

=(aO_l*aO_24M54/eo_max)"054/(Koh.054)4l54/(Ks_ox*l64)4H54 

a(aO_1+aO_24M55/eo_max)4055/(Koh*055)4l55/{Ks_ox+l55)4H55 

=(aO_1+aO_24M56/eo_max)4O56/(K0h+O56)4l56/{Ks_ox+l56)'H56 

=(aO_i»aO_24M57/eo_max)"057/(Koh*057)4IS7/(Ks_ox*l57)4H57 

=(aO_1*aO_24M58/eo_max)4058((Koh*058)4l58/(Ks_ox*l58)4H58 

=(aO_1+aO_24M59/eo_max)4059/(Koh*059)4l59/(Ks_ox+l59)4H59 

=(aO_1+aO_24M80/eo_max)40€0/(Koh+060)4l60/(Ks_ox+l60)4H60 

=(aO_1+aO_24M81/eo_max)"061/(Koh*061)4l61/(Ks_ox+l61)4H61 

=(aO_1+aO_24M62/eo_max)4O62/(K0h+O62)4l62/(Ks_0X*l62)4H62 

=(aO_1+aO_24M63/eo_max)"063/(Koh*063)4l63/{Ks_ox+l63)4H63 

=(aO_1*aO_24M64/eo_max)4064/(Koh*064)"l64/(Ks_ox*l64)4H64 
={aO_1*aO_24M65/eo_max)4065/(Koh*065)4l65/(Ks_ox*l65)‘H65 
=(aO_1*aO_24M66/0O_max)4O66/(Koh+O66)4l66/(Ks_ox4l66)4H66 
=(aCM*aO_2"M67/eo_rriax)4067/(Koh"067)4l67/{Ks_ox»l67)4H67 
=(aO_1*aO_24M68/eo_max)4068/(Koh*068)"l68/(Ks_ox*l68)4H68 
=(aO_1+aO_2"M69/eo_max)4069/(Koh+069)4l69/(Ks_ox+l69)4H69 
s=(aO_HaO_24M70/eo_max)4070/(Koh+070)*l70/(Ks_ox+l70)4H70 
=(aO_1+aO_24M71/eo_max)4071/(Koh+071)"l71/(Ks_ox+l71)4H71 
=(aO_1+aO_24M72/eo_max)4072/(Koh+072)"l72/(Ks_ox+l72)4H72 
=(aO_1+aO_24M73/eo_max)4073/(Koh*073)4l73/(Ks  ox*l73)4H73 
=(aO_1*aO_24M74/eo_max)4074/(Koh*074)4l74/(Ks_ox*l74)4H74 
=(aO_1+aO_24M75/eo_max)4O75/(K0h+O75)4l7S/{Ks_ox+l75)4H75 


={aNO_1+aNO_2'N25/eno_max)'J25/(Kno+J25),l25/(Ks_an*l25)‘H25'Ki/(Ki+025) 
=(aNO_1+aNO_2*N26/eno_max)'J26/(Kno+J26)*l26/(Ks_an+l2B)*H26*Ki/(Ki+026) 
=(aNO_1+aNO_2*N27/eno_max)‘J27/(Kno+J27)*l27/(Ks_an+l27)'H27'Ki/(Kj+027) 
={aNO_1*aNO_2*N28/eno_max)‘J28/(Kno+J28n28/(Ks_an+l28)'H28’Ki/(Ki+028) 
={aNO_1+aNO_2‘N29/eno_max)‘J29/(Kno+J29)M29/(Ks_an*l29)"H29*KI/(Ki+029> 
=(aNO_1+aNO_2’N30/eno_max)M30/(Kno+J30)*!30/(Ks  an*l30)-H30'Kl/(Ki*0301 
=(aN0_1  +aNO_2'N31/eno_max)‘J31/(Kno+J31  p31/(Ks_an*l31  )-H31  •KU(Ki*031 ) 
=(aN0_1  •aNO_2'N32/eno_max)-J32/(Kno*J32)'l32/(KsIan*l32)'H32-Ki/(Ki*032) 
=(aNO_1+aNO_2*N33/eno_max|*J33/(Kno+J33),l33/(Ks_an+l33)'H33‘Ki/(Ki+033) 
-(aNO_1+aNO_2'N34/eno_max)\l34/(Kno+J34)*l34/(Ks_an+l34)'H34'Ki/(KI+034} 
=(aNO_1»aNO_2*N35/eno_max)'J35/(Kno^J35)*l35/(Ks_an*l35)‘H35'KI/(Ki»035) 
=(aNO_1+aNO_2*N36/eno_max}*J36/(Kno*J36)*l36/(Ks_an+l36)*H36,Ki/(KH-036) 
=(aNO_1+aNO_2'N37/eno_max)*J37/{Kno*J37)'l37/{KS_an+l37)‘H37"Ki/(Ki*037) 
=(aNO  1*aNO_2-N38/eno_max)-J38/(Kno*J38)-|38/(Ks  an*l38)’H38'Ki;(Ki*038) 
=(aNO_1.aNO_2-N39/eno_max)-J39/(Kno.J39)-|39/(Ks_an*l39)-H39'Ki/(Ki*039) 
=(aNO_1+aNO_2'N40/eno_max)\l40/(Kno*J40)"W0/(Ks_an+l40)'H40,Ki/(Ki+O40) 
=(aNO_1+aNO_2'N41/eno_max)‘J41/(Kno+J41)‘l41/(Ks_an*l41)*H41’Ki/(Ki«041) 
=(aNO_1taNO_2-N42/eno_max)'J42/(Kno*J42)-l42/(Ks_an*l42),H42'Ki/(Ki*042) 
=(aNO_1*aNO_2-N43/eno_max)'J43/(Kno*J43)-|43/(Ks_an»i43)-H43,Ki/(Ki*043) 
=(aNO_1+aNO_2‘N44/eno_max)'J44/(Kno*J44)'l44/(Ks_an*l44)'H44‘Ki/(Ki+044) 
=(aNO_1*aNO_2‘N45/eno_rnax)'J45/(Kno*-J45)'l4S/(Ks_an*l45)*H45*Ki/(Ki+045) 
=(aNO_1*aNO_2'N46/eno_max)'J46/(Kno*J46)*l46/(Ks_an*l46),H46*Ki/(Ki+046) 
=(aN0_1  *aNO_2*N47'eno_max),J47/(Kno»J47)*l47/(Ks_an*t47)‘H47,Ki/(Ki*047) 
=(aNO_1*aNO_2'N48/eno_max)*J48/(Kno+J48)*l48/(KsIan+148)'H48*Ki/(Ki+048) 
=(aNO_1*aNO_2'N49/eno_max)'J49/(Kno+J49)*l49/(Ks_an+l49)'H49‘Ki/(Ki*049) 
=(aNO_1+aNO_2'N50/eno_max)’J50/(Kno*J50)'l50/(Ks_an+l50)'H50*Ki/{Ki*050) 
=(aN0_1  +aNO_2'N51/eno_max)'J51/(Kno+J51  )'l51/(Ks_an*l51  )'H5rKI/(Ki*051 ) 
=(aNO_HaNO_2'N52/eno_max)’J52/(Kno»J52),l52/(Ks_an+l52)'H52'Ki/(K)*052) 
=(aNO_1*aNO_2,N53/eno_max)'J53/(Kno+J53)'l53/(Ks_an*l53)'H53*Ki/(Kl*053) 
=(aN0_1*aN0  2,N54/eno_max)’J54/(Kno*J54)‘l54/(Ks  an«IS4)-H54'KI/(Ki*054) 
=(aNO_1*aNO_2*N55/eno_max)‘J55/(Kl>o*J55)'l55/(Ks_an»l55)‘H55*Kl/(Kl*055) 
=(aNO_1*aNO_2-N56(eno_max)'J56/(Kno»J56)'l56/(K$_an*l56)'H68'Ki/(Ki‘056) 
=(aNO_1*aNO_2‘N57/eno_max)'J57/(Kno*J57)’l57/(Ks_an*l57)*H57‘Ki/(Ki*057) 
=(aNO_l*aNO_2‘N58/eno_max)'J58/(Kno+J58J*l58/(Ks_an*l58)4H58‘Ki/(Ki+058) 
=(aNO_1+aNO_2*N59/ero_max)‘J59/(Kno+J59)*l59/(Ks_an+l59)"H59‘Ki/(Ki+059) 
=(aNO_l*aNO_2-N60/eno_max)'J60/(Kno*J60)'l60;(Ks_ari*l60)-H60'Ki/(Ki*060| 
=(aN0_1  ♦aNO_2,N61/eno_max)*J61/(Kno+J61  )'l61/(Ks_an*l61  )'H61  'Ki/(Ki*061 ) 
=(aNO_1+aNO_2*N62/eno_max)‘J62/(Kno+J62)'l62/(Ks_an+l62)‘H62*Ki/(Kt+C362) 
=(aNO_1+aNO_2"N63/eno_max)*J63/(Kno+J63),l63/(Ks_an+l63)*H63‘Ky(Ki+063) 
s{3NO_1+aNO_2*N64/eno_max)*J64/(Kno+J64)‘l64/(Ks_an+l64)'H64'K9(Ki+064) 
=(aNO_1+aNO_2,N65/eno_max)*J65/(Kno+J65)‘l65/(Ks_an+l65),H65'Ki/(Ki+065) 
={aNO_1+aNO_2*N66/eno_max)*J66/(Kno+J66)‘l66/(Ks_an+l66)*H66'Ki/(Ki*066} 
=(aNO_1+aNO_2*N67/eno_max)‘J67/(Kno+J67)*l67/(Ks  anH67)-H67,Ki/(KH067) 
=(aNO  1+aNO  2'N68/eno_max)-J68/(Kno*J68|-|68/(Ks  an*l68)'H68'Ki;(Kl*068) 
=(aNO_1+aNO_2*N69/eno_max)*J69/(Kno*J69)"l69/{Ks_an+l69)’H69*Ki/{KH069) 
=(aNO_1+aNO_24N70/eno_max)*J70/(Kno*J70)'l70/(Ks_an+l70)*H70'Ki/(KH070) 
=(aNO_1+aNO_2,N71/eno_max)‘J71/(Kno*J71)*l71/(Ks_an+l71!‘H71'Kt/(Kt+071) 
=(aNO_1*aNO_2*N72/eno_max)‘J72/(Kno*J72)’172/(Ks_an+l72)‘H724Kl/(10+072) 

=(aNO_1*aNO_2*N73/eno_max)'J73/(Kno*J73)-|73/(Ks_an*l73)'H73‘Ki/(K(*073) 

=(aNO_1+aNO_2*N74/eno_max)'J74/(Kno+J74)"l74/(Ks_an*l74)4H74*KI/{Ki+074) 

=(aNO_1+aNO_2*N75/eno_max)‘J75/(Kno+J75)*l75/(Ks_an*175)‘H75‘Ki/(Ki*075) 


_ox)‘O25/(K0h+O25)*l25/(Ks_ 
ox)*026/(Koh+026)*l26/(Ks_ 
0X)*O27/(K0h+O27)*l27/(Ks 
_ox)*028/(Koh+028)*l28/(Ks_ 
0x)'O29/(K0h+O29)*l29/(Ks_ 
ox)*030/(Koh+030)*l30/(Ks_ 
ox)*031/(Koh+031  )‘I31/(KS_ 
_ox)'032/(Koh+032)*l32/(Ks_ 
'n_ox)*033/(Koh+033)‘l33/(Ks_ 
'n_ox)'034/(Koh+034)*i34/(Ks_ 
'n  ox)‘035/(Koh+035)‘l35/(Ks 
’n_0X)'O36/(Koh*O36)'l36/(Ks_ 
'n~0X)'O37/(K0h*O37)‘l37/(Ks_ 
'n_ox)'038/(Koh+038)'l38/(Ks_ 
'n  ox)'039/(Kon*039)‘l39/(Ks 
vn_ox)*040/(Koh*040)*l40/(Ks_ 
'n_ox)*041/(Koh*041  )-|41/(Ks_ 
'n_ox)’042/(Koh+042)'l42/{Ks_ 
'n_ox)*043/(Koh+043)’l43/(Ks_ 
’n_ox)’044/(Koh»044)*l44/(Ks_ 
'n~ox)’045/(Koh*045)*t45/(Ks_ 
'n_ox)‘046/(Koh+046)’l46/(Ks_ 

'n_ox)-O47/(K0h*O47)’l47/(Ks_ 
'n_ox)‘O48/(K0h*O48)*l48/(Ks_ 
'n_ox)'049/(Koh*049)'l49y(Ks_ 
'n_ox)‘050/(Koh*050)'l5Q/(Ks_ 
'n  ox)-051/(Koh*051)'l51/(Ks 
'n_ox)-052/(Koh‘052)'l52/(Ks_ 
'n_ox)‘053/(Koh*053)*l53/(Ks_ 
'n_ox|*054/(Koti*054(*l54/(Ks 
■n  ox|*055/(Koh*055)M55/(Ks 


'n_ox)‘O57/{K0h+O57)*l57/(Ks_ox+l57)*H57 
'n_ox)*058/(KcMl*058)'l58/(Ks_ox*158}*H58 
'n  ox)*059/(Koh*059)M59/(Ks_ox*l59)'H59 
'n_ox)'060/(Koh*060)'l60/(Ks_ox*l60)'H60 
"n_0X)‘O61/(Koh*O61)‘l61/(Ks_QX*l61)’H61 
'n_ox)*062/(Koh+062),l62/(Ks_ox+l62)*H62 
■n  ox)'063/(Koh*063)*163/(Ks_ox‘l63)'H63 
>n  ox)'064/(Koh+064)*l64/(Ks  0X*I64|'H64 


(+I25)*H2S 

(*I26)‘H26 

0X+I28)*H28 
(+I29)*H29 
C+I30)*H30 
_0X+I31  )*H31 

’ox*I34)'H34 

’ <*I36)‘H36 
C+I37)*H37 
0X*I38)‘H38 

OX+I40)"H40 

,0X+I41)*H41 

<*142)’H42 

0X*I44)'H44 

"ox»I48)*H48 

0X»I52)'H52 


ox)'066/(Koh*066)'l66/(Ks_ 
,OX)’067/(Koh+067)‘l67/(Ks_ 
ox)’068/(Koh*068)-l68/(Ks_ 
oxf069/(Koh*069)-l69/(Ks 
ox)*07Q/(Koh+070)‘l70/(Ks 
px)‘071/(Koh*071  )*I71/(Ks_ 
px)'072/(Koh*072)‘l72/(Ks_ 
px)-073/(Koh*073)-|73/(Ks_ 
_OX)'074/(Koh*074)*l74/(Ks^ 
ox)*075f(Koh«075)*l75/(Ks 


_0X+I68)‘H68 

OX*I70)‘H70 

0X*I71)'H71 

_0X*I72)‘H72 

_0X+I73)*H73 


=(p+(N25/eno_.max)An_anoxAmumax_an)*Ki/(Ki+025)"J25/{Kno+J25)'l25/(Ks_an+l25)'H25 

=(pA(N26/eno_max)An_anox*mumax_an)'Ki/(Ki*026)\l26/{Kno+J26)'l26/(Ks_an+l26)'H26 

=(p+(N27/eno_max)An_anox*mumax_an)'Ki/(Kt+027)‘J27/{Kno+J27)'l27/(Ks~an+l27)‘H27 

=(p*(N28/eno_max>An_anox'mumax_an)*Ki/(Ki+028)*J28/(Kno+J28)'l28/(Ks_an+l28)*H28 

=(p+(N29/eno_max)An_anox"mumax_an)'Ki/(Ki+029)\l29/(Kno+J29)'l29/(Ks_ar>+l29)*H29 

=(p+(N30/eno_max)An_anox*mumax_an)*Ki/{Ki+030)*J30/(Kno+J30)‘l30/(Ks_arvH30)’H30 

=(p*(N31/eno_max)An_anoxAmumax_an)-Ki/(Ki*031)AJ31/(Kno*J31)Al31/(Ks_an*131)AH31 

=(p*(N32/eno_max)An_anoxAmumax_an)'Ki/(K.*032)'J32/(Kno»J32)'l32/(KsIan-H32)-H32 

=(p*{N33/eno_max)An_anox‘mumax_an)'Ki/{Ki+033)*J33/(Kno*J33)'l33/(Ks_an+l33)'H33 

=(p+(N34/eno_max)An_anox'mumax_an)‘Kr/{Kt*034)‘J34/(Kno+J34)*l34/(Ks_aivH34)’H34 
=(p*(N35/eno_max)An_anox-mumax_an)-Ki/(Ki*035)'J35/(Kno*J35)-|35/(Ks_an*l35rH35 
=(p+{N36/eno_max)An_anox‘mumax_an)AKi/(Ki+036)\J36/(Kno+J36)*l36/(Ks_an+l36),H36 
=(p+(N37/eno_max)An_anox’mumax_an)*Ki/(Ki+037)\J37/(Kno*J37)‘l37/(Ks  an*l37)AH37 
=(p+(N38/eno_max)An_anox*mumax_an>*Ki/(Ki*038)*J38/(Kno+J38)*l38/{Ks_an+l38)*H38 
-(p+(N39/eno_max)An_anox*mumax_an)'Ki/(Ki+039)’J39/(Kno+J39p39/(Ks__an*l39)AH39 
=(p+(N40/eno_max)Ar_anoxAmumax_an)*Ki/(Ki+040)'J40/(Kn&»J40)*l40/(Ks_an+l40)*H40 

=(p+{N4Veno_max)An_anoxAmumax_an)AKi/{KI-*041)’J41/(Kno*J41)‘l41/{Ks_an+l41)‘H41 

=(p+(N42/eno_max)An_anox*mumax_an)*Ki/(KI+042)'J42/(Kno+J42)'l42/(Ks_an+l42)*H42 

=(p+(N43/eno_max)An_anox'mumax_an)*Ki/(KI+043)‘J43/(Kno+J43)*l43/(KsIan+l43|,H43 

=(p+(N44feno_max)An_anox*mumax_an)*Ki/(Ki*044)‘J44/(Kno+J44)’l44/(Ks_an+l44)‘H44 

=(p*(N45feno_max)An_anoxAmumax_an)AKi/(Ki*045)'J45/(Kno*J45)‘l45f(Ks_an*l45)AH45 

=(p*(N46/eno_max)An_ancix"mumax_an)-KI/(Ki-046)'J46/(Kno*J46)'l46/(Ks_ari*l46)'H46 

=lp+(N47/eno_max)An_anox‘mumax_an)*Ki/(Ki-*047)‘J47/(Kn&AJ47)Al47/(Ks_an*l47)*H47 

a(p+(N48/eno_max)An_anox*mumax_an)‘Ki/(Ki+048)\J48/{Kno+J48)*l48/(Ks_an*l48)'H48 

={p*(N49/eno_max)An_anox'mumax_an)'KI/(KI*049)-J497(Kno*J49)-l49/(Ks_an*l49|-H49 

-(p*(N50/eno_max)An_anox*mumax_an|'Ki/(Ki»050)*J50/(Kno*J50)‘l50/(Ks_an*l50)*H50 

={p+(N51/eno_max)An_anox'mumax_an)*Ki/(Ki+051),J51/(Kno+J51)*l5V(Ks_an+151)*H51 

=(p*(N52/eno_max)An_anox'mumax_an)'Ki/(Ki*052)-J52/(Kno*J52)Al52/(Ks  an*l52)‘H52 
=(p*(N53/eno_max)"n_anox‘mumax_an|‘Ki/(Ki*053)*J53/(Kno*J53)*l53/(Ks_an*l53}’H53 
=(p*(N54/eno_max)An_anox'mumax_an),Ki/(KH054)'J54/(Kno»J54)Al54;(Ks_aii»l64)AH54 
={p+(N55/eno_max)An_anoxAmumax_an)'Ki/(Ki+055)*J55/(Kno+J55)*l55/(Ks_an+155),H55 
a{p+(N56/eno_max)An_anox'mumax_an)"Ki/(Ki-*056)'J56/(Kno+J56)*l56/{Ks  an*l56|-H56 
a(p+(N57/eno_max)An_anox'mumax_an)*Ki/(Ki+057)\l57/(Kno+J57)*l57/(Ks_an+l57)*H57 
-(p*(N58/eno_max)An_anox'mumax_an}*Kl/(Ki+058)'J58y(Kno*J58)'l58/(K$Ian*l58)-H58 
=(p‘(N59/eno_max)An_anox'mgmax_anrKi((Ki*059)'J59/(Kno‘J59)-l59/(Ks_an*l59)-H59 
a(p*(N60/eno_max)An_anox’mumax_an)'Ki/(Ki+060>"J60/(Kno+J60)*l60/(Ks_an+l60)'H60 
a(p*(N61/eno_max)An_anox‘mumax_an),Ktf(Ki+061)'J61/(Kno+J61)'l61/(Ks_an+l61)'H61 
=(p*(N62/eno_max)An_anox-mumax_an)'Kl/(Ki*062)-J62((Kno*J62)Al62/(Ks_an*l62)'H62 
=|p*(N63/eno_max)An_anoxAmumax_an)*Kt/(Ki+063)\J63/(Kno+J63)'l63/(Ks_an+l63)'H63 

=(p*(N64/eno  max)An  anox'mumax_an)'Ki/(Ki-t064)-J64/(Kno*J64)'l64/(Ks  an*l64)'H64 

=(p*(N65/eno_max)An_anox‘mumax_an)*Ki/(Ki+065)*J65/(Kno+J65)Al65/(Ks_an+l65)"H65 

=(P*(N66/eno_max)An_anox'mumax_an)'Ki/(Ki+066)*J66/(Kno+J66n66/(Ks_an+l66)'H66 

=(p+(N67/eno_max)An_anox'mumax_an)‘Ki/(Kt+067)‘J67/(Kno+J67)Al67/(Ks_an+l67)*H67 

=(p+(N68/eno_max)An_anox'mumax_an)*Ki/(Ki+068)’J68/(Kno+J68n68/(Ks_an+l68)'H68 

=(p*(N69/eno_max)An_anox*mumax_an)*Ki/(Ki*069|*J69/(Kno*J69)*l69/(Ks_an*l69)*H69 

=(p*{N70/eno_max)An_anox*mumax_an)'Ki'(Ki+070)AJ70/(Kno+J70)*l70/(Ks_an+l70)'H70 

=(p+(N71/eno_max)An_anox'mumax_an)*Ki/(Ki+071)*J71/(Kno+J71)'l71/(Ks_an+l71)*H71 

a(p+(N72/eno_max)An_anox'mumax_an)*Ki/(Ki*072)M72/(Kno+J72)'l72/(Ks_an+l72)AH72 

=(p+(N73/eno_max)An_anox,mumax_an),Ki/{Ki*073)'J73/(Kno*J73)'l73/(Ks_aft+l73)-H73 

=(p+(N74/eno_max)An_anox'mumax_an)*Ki/(Ki+074)'J74/(Kno*J74)*l74/(Ks_an+l74)‘H74 

a(p+(N75/eno_max)An_anox*mumax_an)‘KI/(Ki+075)\l75/(Kno+J75)'l75/(Ks_an+l75)*H75 
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=R25+S25-(b+D)"H2! 

=R26*S26-(b*DfH2l 

=R27+S27-(b+D)*H2' 

=R28+S28-{b+D)*H2l 

=R29+S29-{b+D)*H2! 


=R31*S31-(b*D)*H3 

=R33+S33-(b+D)*H3! 

*'34+S34-{b+D)*H3‘ 


=R37+S37-{b*D)*H37 

=R39*S39-(b*D)*H39 

=R40*S40-(b+D)*H40 

=R41*S41-(btD)*H41 

=R42+S42-(b*D)*H42 

=R43»S43-<b*D)*H43 

=R44+S44-(b*D)"H44 

=R45+S45-(b+D)*H45 

=R46+S46-(b+D)*H46 

=R47+S47-(b+D)"H47 

=R48+S48-(b+D)*H48 

=R49*S49-(b*D)*H49 

=R50+S50-(b+D)*H50 

=R51»S51-(b*D)*H51 

=R52*S52-<b*D)*H52 

=R53*S63-(b*D)*H53 


=(*1/Yc_ox)*R25-(1/Y c_an)*S25+(1 000-I25)*D 
=(-1/Yc  ox)*R26-(1/Yc_an)*S26+(1000-l26)*D 
=(-1/Yc_ox)‘R27-(l/Yc_an)*S27-*(1000-l27|*D 
-(-17r'c_oxfR28-(l/Yc_an)*S28+f 1000-1281*0 
=(-1/Yc_ox)*R29-(1/Yc_an)*S29*(1000-l29)*D 
=(-1/Yc_ox)‘R30-(1/Yc_an)‘S30*(1000-l30)*D 

=(-1/Yc_ox)*R31-(1/Yc_an)*S31+(1000-l31)*D 
=(-1/Yc  ox)*R32-(1/Yc_an)*S32+(  1 000-I32)‘D 
=(-1/Yc_ox)*R33-(1/Ye_an)*S33*(1000-l33)*D 
=(-1/Yc_ox)*R34-(1/Yc  an)*S34*(1000-l34)*D 
=(-1/Yc_ox)*R35-(1/Yc_an)*S35+(1000-l35)*D 
=(-1/Yc_OX)*R36-(1/Yc  an)*S36+(1000-136)*D 
=(-1/Yc_ox)*R37-{1/Yc_an)*S37*<1000-l37)’D 
=(-1  /Y c_ox)*R38-(  1 /Y c_an)*S38+( 1 000-I38)*D 
~{-17Yc_OX)*R39-(l/Yc_anl*S39*(1000-l391*D 
=(-Wc_ox)*R40-<1/Yc_ari)*S40*(1000-l40)*D 
=(-lA'c_ox)*R41-{lA'c_an)‘S41-*(  1000-141 1*0 
=(-1/Yc_ox)*R42-<1/Yc_an)*S42*(1000-l42)*D 
=(-17Yc_ox)*R43-(1/Yc_an)*S43*(1000-l43)*D 
=(-1  /Y c_ox)*R44-(1  /Y 0_an)*S44+( 1 000-M4)*D 
=(-177c_OX)*R45-(17Yc_an)*S45*(1000-l45)*D 
=(-l/Yc_oxl*R46-(l/Yc_an)‘S46*(  1000-1461*0 
=(-1/Yc_ox|*R47-(UYc_an|*S47*(1000-l47)*0 
=(-1/Yc_ox)‘R48-(1/Yc_an)*S48+(1000-l48)*0 
=(-1/Yc_ox)*R49-(11Yc_an)*S49*(1000-l49)*D 
=(-1/Yc_ox)*R50-(1/Yc_an)*S50+(1000-l60)*D 
=(-1/Yc_ox)*R51-(17Yc_an)*S51*(1000-IS1)*D 
=(-1/Yc_ox)*R52-(1/Yc__an)*S52*{1000-l52)*D 
=(-1/Yc_ox)*R53-(1/Yc_an)*S53+{1000-l53rD 
=(-1/Yc  ox)*R54-(1/Yc_an)*S54+(1000-l54)*D 
=K-1/Yc_OX)*R55-(1/Yc_an)*S55»(100O-l55)*D 
=(-17Yc_ox)*R56-(17Yc_an)*S56+(1000-l66)*D 
=(-1/Yc  ox)*R57-(1/Yc_an)*S57*(1000-l57)‘D 
=(-1/Yc_OX)*R68-(11Yc_an)*S58*(1000-l58)-D 
=(-1/Yc_ox)*R59-<1/Yc_an)*S59+(1000-l59)*D 
B(-1A'c_ox)*R60-(1/Yc_an)*S60+(1000-l60)*D 
=(-1/Yc_OX)*R61-(1/Yc_an)*S61*(1000-l61)*D 
={-1  /Yc_ox)*R62-<  1 /Yc_an(*S62+(1 000462  )*D 
=(-lA'c_ox)*R63-(l/Yc_an)*S63.(1000-l63)*D 
~(-1/Y  C_OX)*R64-{  1 /Y c_anl*S64+( 1 000-164 )*D 
-(-17Y c_ox)*R65-(  1 /Y c_an)*S65+( 1 000-I65)*D 
=(-1/Yc_ox)*R66-{lA'c_an)*S66*(1000-l66)*0 
=(-11Yc_ox)*R67-(1/Yc  an)*S67+(1000-I67|*D 
=(-1/Yc_OX)*R68-(1/Yc_an)*S68*(1000-l68)'D 
=(-17Yc_OX)*R69-(1/Yc_an)*S69+(1000-l69)*D 
=(-1/Yc_ox)*R70-(1/Yc_an)*S70«(1000-l70)*D 
=(-1/Yc_ox)*R71-(1/Yc_an|*S71+(1000-l71)*0 

=(-1/Yc  Ox)*R72-(17YC_an)*S72»(1000-l72)*0 
=(-1/Yc_OX)*R73-(1/Yc_an)*S73*(1000-l73)*D 
=(-1/Yc_ox)*R74-(1/Yc_an)*S74+(1000-l74)*D 
=(-1/Yc_ox)*R75-(l/Yc_an)*S75»(1000-l75)*0 
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=-((Wc_an)/(2.86-Yc_an)|-S25 
=-{(1  A'c_an)/(2  86'Yc  an))'S26 
H(1/Yc_an)/(2.86-Yc_an)|-S27 
=-((1'Yc_an)/(2  86-Yc_an))‘S28 
=-((1/Yc_an)/(2.86*Yc_an)}'S29 
=-((1/Yc_an)((2.86’Yc_an))*S30 
=-((1/Yc_an)/(2  86'Yc_an))‘S31 
=-((1'Yc_an)/(2  66'Yc_an))-S32 
=-((1/Yc_an)/(2  86'Yc_an))-S33 
=-((1/Yc_an)/(2.86'Ycj>n))‘S34 
=-«1/Yc_an)/(2  86-Yc_an))‘S35 
=_((1/Yc_anJ/(2.86*Yc_an))*S36 
=-((1/Yc_an)((2  86-Yc_an)rS37 
=-((1/Yc_an)/{2  86"Yc_an))*S3B 
H(1/Yc_anV(2.86-Yc_an))1S39 
=mWc_anV(2.86-Yc_an))-S40 
H(1'Vc_anW2.86-Yc_an))-S41 
=-((1/Yc_an)/(2.86*Yc_an))'S42 
=K(1A'c_anV(2.86'Yo_an))-S43 
H(17Yc  an|/(2  86-Yc_an))'S44 
=-f(177c_anK(2  86*Yc_an|)‘S45 
=K(1/Yc_anV(2.86-Yc_an))-S46 
H(Wc_anV(2.86-Yc_an))-S47 
H(Wc_anV(2.86-Yc_an))-S48 
=-({1/Yc_any(2.86‘Yc_an))‘S49 
--<(1770  an)/(2.86'Yc  an))‘S50 
H(1'Yc_an)/(2B8-Yc_an))-S51 
=-«1/Yc_an)/(2.86*Yc_an))*S52 
=-«  1 TYc  an)/(2.86‘Yc_an))*S53 
=-((17Yc_an)/(2.86*Yc_an))‘S54 
=-((17Yc_an)/(2.86*Yc_an))*S55 
H(1(Vc  an)/(2.86'Yc.  an))‘S56 
=-((17Yc_an)/[2.86*Yc_an))*S57 

=-((1(Yc_an)/(2.86*Yc_an)),S58 
=-((Wc_an)/(2.86-Yc_an))-S59 
=-((1/Yc_an)/(2  86*Yc_an))‘S60 
=-((1(Yc_an)/(2.86*Yc_an))'S61 
=-((17Yc_an)/(2.86*Yc_an))*S62 
--((17Yc_an)/(2  86*Yc_an))‘S63 
=-((1/Yc_an)/(2.86*Yc_an))*S64 
=-(|1/Yc_an)/{2.86‘Yc_an))*S65 
=-((1'Yc_an)/{2.86*Yc_an))*S66 
=-«1A'c_an)/(2.86-Yc_an))-S67 
=-((1'Yc_an«2.86-Yc_an))-S68 
=.((WYc_anV(2.86-Yc_an))-S69 
=-((  1/Yc_an|/(2.86-Yc_an))'S70 
=-((Wc_anV(2.86-Yc_an»'S71 

=-!(  1/Yc_any(2.86*Yc_anJ)'S72 

H<1/Yc_anV(2.86-Yc_an))-S73 

=-<(1"cIan)/(2.86-YcIan))-S75 


=P25-(bO+D)‘K25 

=P26-(bO+D)*K26 

=P27-(bO+D)'K27 

=P28-(bO+D|*K28 

=P29-(bO+D)'K29 

=P30-(bO+D)*K30 

=P31-(bO*D)*K31 

=P32-(bO+D)*K32 

=P33-(bO+D)*K33 

=P34-(bO+D)‘K34 

=P35-(bO+D)*K35 

=P36-(bO+D)'K36 

=P37-(bO+D)'K37 

=P38-(bO+D)'K38 

=P39-(bO+D)*K39 

=P40-(bO+D)'K40 

=P41-<bO+D)*K41 

=P42-(bO+D)*K42 

=P43-<bO+D)‘K43 

=P44-{bO+D)*K44 

=P45-<bO+D)*K45 

=P46-<bO*D)'K46 

=P47-(bO+D)*K47 

=P48{bO*D)'K48 

=P49-(bO+D)’K49 

=P5<HbO*D)'K50 

=P51-(bO+D)*K51 

=P52-(bO+0)‘K52 

=P53-<bO+D)*K53 

=P54-(bO*D)‘K54 

■pe&Kbo+oj'KM 

=P56-(bO+D)*K56 

=P57-(bO*D)'K57 

=P59-(bO*D)'K59 

=P60-(bO*D)'K60 

=P61-(bO*D)‘K61 

=P62-(bO*D>*K62 

=P63-(bO+D)’K63 

=P64-(bO+D)'K64 

=P66-(bO+D)*K66 

=P87-|bO*D)-K67 

=P68-(bO*D)'K68 

=P69-<bO+D)'K69 

=P70-(bO+D)*K70 

=P71-<bO+D)*K71 

aP72-<bO+D)*K72 

=P73-(bO+D)'K73 

=P74-(bO+D)‘K74 

=P75-<bO*D)"K75 


=Q25-<bNO+D)*l25 

=026-<bN0*DCl26 

=Q27-<bNO*D)-L27 

=Q28-{bNO+D)'L28 

=029-<bN0+D)*t29 

=O30-(bNO+D)*L30 

=Q31-<bNO+D)*L31 

=Q32-<bNO+D)*L32 

=Q33-(bNO+D)'L33 

=Q34-<bNO*D)'L34 

=035-(bN0*D)'L35 

=Q36-(bNO+D)‘L36 

=Q37-(bNO+D)*L37 

=Q38-(bNO*D),L38 

=Q39-(bNO+D)‘L39 

=Q40-(bNO*D)*L40 

=Q41-(bNO*D)*L41 

=042-<bN0»D)'L42 

aQ43-(bNO*D)'L43 

=Q44-(bNO*D)‘L44 

=Q46-(bNO+D)'L45 

=Q46-(bN0*0)"L46 

=Q47-(bNO*D)-L47 

=Q48-(bNO*D)-L48 

=Q49-(bNO*D)‘L49 

=Q50-(bNO*D)'L50 

=Q51-(bNO*D)-L51 

=Q52-(bNO+D)‘L52 

=Q53-(bNO*D)-L63 

=Q54-(bNO+D)*L54 

=Q55-(bNO»D)-L55 

=Q56-(bNO*D)*L56 

=Q57-(bNO+D)*L57 

=Q58-<bNO+D)*L58 

=059-<bN0*D)‘L59 

=Q60-<bNO+D)'L60 

=Q61-<bNO*D)*L61 

=Q62-(bNO+D)‘L62 

aQ63-(bNO+D)'L63 

=Q64-(bNO+D)'L64 

=Q65-(bNO+D)'L65 

=Q68-(bNO+D)'L66 

=Q67-<bNO+D)'L67 

=Q68-(bNO*D)'L68 

=Q69-(bNOD)’L69 

=Q70-(bNO*D)*L70 

=Q71-(bNO+D)*L71 

=Q72-(bNO+D)*L72 

=Q73-(bNO*D)‘L73 

=Q74-(bNO+D)*L74 

=Q7S-(bNO*D),L75 


APPENDIX  B 

SEMI-LOGARITHMIC  PLOTS 


The  semi-logarithmic  plots  of  data  from  the  anoxic  phase  of  experiments 
discussed  in  Chapter  5 are  shown  in  this  section.  Growth  of  cultures  using  acetate  or 
malate  as  carbon  substrate  resulted  in  two  straight  lines  of  different  slopes  (which 
represent  two  different  growth  rates).  Growth  of  Faracocciis  panloirophus  using 
caproate  resulted  in  one  straight  line  with  one  slope  throughout  the  entire  anoxic 
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APPENDIX  C 

ADDITIONAL  MODEL  FITS 


biomass  mg/LCOD 
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